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1.  INTRODUCTION 

For  the  convenience  of  analysis,  mooring  line  tension  under  the  excitation 
of  wind,  wave  and  current  may  be  divided  into  a  steady  state  component  and  a 
dynamic  component.  The  steady  state  mooring  line  tension  is  defined  as  the 
force  induced  by  the  mooring  line  geometry,  gravity,  average  current  and  aver¬ 
age  wind,  and  the  dynamic  tension  by  wave  and  wind  gusts.  The  steady  state 
analysis  has  been  treated  extensively  and  the  solution  can  be  achieved  in  a 
relatively  small  amount  of  computer  time.  The  dynamic  analysis  has  not  been 

developed  to  the  point  where  a  lew  cost  computer  program  with  good  accuracy 

(2)  (3) 

is  available.  The  works  of  Paquette  and  Henderson,  Wilson  and  Garbaccio, 
Reid/4^  Kaplan  and  Ralf,^  Nath/6^  and  Brainard^  are  typical  of  the 
studies  of  the  dynamic  response  of  single  point  moorings  made  in  recent  years. 

Intensive  reviews  of  the  literature  on  the  response  of  various  cable 
systems  under  hydrodynamic  loading  are  presented  by  Casarella  and  Parsons. ^ 

A  time  domain  analysis  is  most  desirable  for  the  solution  of  mooring  line 
tension  under  random  waves.  The  solution  technique  may  be  divided  into  two 
types:  the  digital  computer  approach  through  tie  use  of  the  method  of  charac¬ 
teristics  and  the  analog  computer  approach.  The  latter  has  been  discussed  in 
Kaplan  and  Raff.^  The  basic  formulation  of  the  field  equations  based  on  the 

method  o^  characteristics  was  presented  by  Reid^4^  and  Nath,^^  and  a  digital 

(  0  ) 

computer  program  was  developed  by  Nath. 

(9) 

By  supplying  a  simulated  random  wave  as  the  excitation  to  the  buoy 
system,  the  random  stress  history  at  any  point  of  the  mooring  line  can  be 
obtained.  Based  on  the  simulated  wave  and  the  stress  history,  the  sutocovari- 


ance  function,  the  cross  covariance  function,  stress  peak  distribution  curve. 


and  the  average  frequency  can  be  derived.  However,  each  of  the  steps,  i.e. 
wave  simulation,  solution  program  and  statistical  analysis,  involves  extensive 


r^'^^'V’-  -r^wv- 


computer  time.  The  high  cost  and  special  knowledge  required  by  this  method  may 
make  it  too  expensive  for  normal  design  use. 

In  the  light  of  the  difficulties  encountered  in  a  complete  time  domain 
analysis,  it  is  necessary  to  search  for  an  analytical  solution  in  the  frequency 
domain  in  order  that  the  dynamic  response  of  a  buoy  system  under  random  input 
can  be  performed  in  a  small  amount  of  computer  time.  Additionally,  little 
knowledge  of  time-series  analysis  is  needed.  However,  the  frequency  domain 
analysis  is  only  applicable  to  a  linear  structural  system, and  consequently  the 
non  -linear  system  has  to  be  linearized. 

The  behavior  of  the  mooring  system  is  heavily  dependent  on  the  hydro- 
dynamic  drags  on  the  surface  buoy  and  the  mooring  rope.  As  will  be  shown  later, 
the  tangential  mooring  line  hydrodynamic  drag  is  a  major  damping  factor  in  the 
response  of  a  deep  sea  mooring  line  subject  to  oscillatory  longitudinal  motion 
at  one  end.  However,  neither  theoretical  solutions  nor  experimental  data  are 
available  for  the  estimation  of  the  tangential  drag  coefficient  of  a  rope  under 
oscillating  motion. 

A  review  of  the  drag  coefficients  on  the  mooring  rope  and  the  surface 
buoy  is  presented  in  Appendix  A. 


2.  STEADY  STATE  MOORING  LINE  TENSION 

In  two  dimensional  analysis ,  the  steady  state  tension  in  the  mooring  line 
is  defined  as  the  force  induced  by  geometry,  gravity,  and  the  c.oplanar  average 
current  and  wind.  This  is  considered  the  best  approximation  of  the  mean  stress 
in  the  mooring  line  subject  to  ei  definite  combination  of  wind,  wave  and  current, 
and  the  coplanar  assumption  gives  a  conservative  solution. 

Consider  the  two  dimensional  free  body  of  an  elemental  length  of  the 
mooring  line  as  shown  in  Fig.  1. 
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Figure  1.  Free  Body  of  a  Mooring  Line  Element 

The  equilibrium  equations  in  normal  and  tangential  directions,  after  neglecting 
the  second  order  terms,  are 

(Fdn  +  W  sin<i>)  AS  =  T  A<}>  (1) 

*nd  (Fdt  -  W  cos<j>)  AS  =  AT  (2) 

where  W  is  the  weight  of  the  rope  in  water  per  unit  length, 

Fp^  and  F^  are  as  in  Equation  (A. 2)  and  (A. 3). 


Treating  the  mooring  line  as  a  series  of  finite  chords,  the  solution  of 
(1)  and  (2)  can  be  approximated  by  incremental  numerical  integration.  The 
computer  program  is  presented  in  Appendix  C. 

The  program  can  handle  the  compound  mooring  line,  made  of  wire  rope  and 
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3.  DYNAMIC  MOORING  LINE  TENSION 

3.1  Model  Simplifications 

For  reasons  stated  in  the  Introduction, the  dynamic  mooring  line  tension 
:.'  -blem  is  solved  using  the  frequency  domain  approach.  Since  the  frequency 
domain  analysis  is  only  applicable  to  a  linear  structural  system,  the  buoy 
system  has  to  be  simplified. 

The  basic  assumptions  of  the  buoy  model  considered  in  this  x*ork  are  as 
follows : 

(1)  The  buoy  is  a  surface  follower  buoy  so  that  the  buoy  response  spec¬ 
trum  can  be  considered  to  be  the  same  as  the  wave  spectrum. 

(2)  The  mooring  line  is  taut  and  can  be  treated  as  a  straight  string. 

The  dynamic  force  in  the  mooring  line  due  to  the  horizontal  movement  of  the 
surface  buoy  under  the  action  of  wave  and  fluctuating  wind  is  negligible  com¬ 
pared  to  that  due  to  the  vert'  il  motion  of  the  buoy. 

(3)  The  tangential  hydrodynamic  drag  on  the  mooring  line,  which  is  pro¬ 
portional  to  velocity  squared,  can  be  linearized  through  a  principle  of 
equivalent  linearization. 

(4)  The  internal  damping  of  the  mooring  ropes  is  linear  and  the  dynamic 
stress  strain  relation  under  sinusoidal  motion  is  given  by 

o  =  (E  +  iE„)e 

1  i. 

Depending  on  the  material  behavior  under  dynamic  loading,  the  linear  damping 
material  may  be  represented  by  several  types  of  mathematical  models.  One  of 
the  models  listed  in  Table  1  may  be  used  to  describe  the  behavior  of  mooring 
lines  and  the  models  are  incorporated  into  the  solution  program  DYNSIN  and 
DYNRAN  as  described  in  the  Appendices. 

(5)  Stress  due  to  strumming  is  neglected. 
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be  expected.  However,  the  assumption  will  provide  an  upper  bound  solution  to 
the  dynamic  tension  of  the  mooring  line  provided  the  natural  frequency  of 
the  heave  motion  of  the  buoy  is  far  from  the  effective  wave  frequency.  In  the 
case  of  a  stationary  buoy,  e.g.  spur  buoy,  the  dynamic  tension  may  be  con¬ 
sidered  as  negligible.  The  second,  fourth  and  fifth  assumptions  introduce 
negligible  error,  as  discussed  in  Appendix  B.  The  error  due  to  the  lineariza¬ 
tion  of  the  hydrodynamic  drag  will  depend  on  the  degree  of  non-linearity.  The 
distortion  of  the  result  may  be  negligible  at  a  low  sea  state  and  significant 
at  a  high  sea  state. 

No  error  bound  is  available,  and  the  accuracy  can  only  be  checked  by 
experimental  data  as  discussed  in  Appendix  B. 

3.2  Dynamic  Mooring  Line  Tension  Under  Sine  Wave 

A  deep  sea  mooring  line  with  instrument  packages  inserted  in  it  is  ideal¬ 
ized  as  Fig.  2. 

The  equilibrium  condition  of  the  free  body,  after  neglecting  the  second 
order  terms,  leads  to: 
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Fig.  2-  Idealized  Mooring  Line  and  Free  Body  Diagram  of  the  Rope  Element 


Introducing  the  relation  a-  (E^  +  iE2)e  =  (E^^  +  i£2>  ~ 

into  Equation  (3)  and  dividing  by  AE^,  then 

32u  .e2  92u  CDTirDpW  3u  3u  pr  32u  _ 

3x2  E1  3x2  "  2AE1  3t  ‘  3t  "  E1  3t2  " 

Equation  (4)  may  be  linearized  to 

■  i_  32u  p  3u  1  32u  _ 

3x2  *  3x2  DR  9t  a2  3t2  ' 

where 

_h  ,  _h 

T  “  g  »  a  **  ~  is  the  equivalent  linear  damping  parameter 

1  ^r 

derived  in  (B.27). 

Separating  variables  in  the  form 

U(x,t)  =  X(x).eiut 

where  X(x)  is  complex,  then  the  spatial  part  of  Equation  (5)  becomes 


ti  -5',^!;  -!ij£§&4k£ 


(1  *  iOX,  x  +  (— —  -  iuC-._)X  =  o 
a2  0R 


Thin  has  the  solution 


X(x)  =  (R^  +  i^)  sin  (at  iR)x  +  (R2  +  il2)  cos  (a  +  i$)x  (8) 

where  R  ,  1^,  R2,  I2  are  real  constants  to  be  determined  from  boundary  and  con¬ 
tinuity  conditions,  and 

T [<l)2-  ,CDR“]  rtC(f>2-  ’CDR“]2  +  CCDR“«<I>^2^  T 

a=  -  ( 9 ) 

2(1  +T2) 

-t(l»2-  ^DR”3  +{C(f>2-  ^DR^  *  tCDR“  *  '<&*■?*' f 

p —  -  V.JLO; 

L  2(1  +T2)  J 

The  values  of  a,  rand  will  depend  on  the  material  model  and  the  mater¬ 
ial  constants.  They  are  constants  for  the  two  parameter  models  and  frequency 
dependent  functions  for  the  three  parameter  rodeis. 

The  relation  between  L , ,  E2  and  the  material  constants  for  various  models 
are  shown  in  Table  1. 

The  displacement  is  now 

U(x,t)  =  [(R^  +  il^)  sin  (a  +  i8)x  +  (R2  +  il^)  ccs  (a-  iB)xje^tl)t  (11) 

The  strain  is 

r  ^R-j+ilj^fo+iiOcosGx+ifOx-  (R2+^2^a+^^^sill^a+i®^y'^eit0t  (12) 

and  the  stress  is 

"x  "  <E1  +  iE2>  i  U3) 

Considering  only  the  real  parts  and  rearranging  gives  the  displacement 

U(x,t)  =  A^cos  wt  +  A2  sin  wt  =  U^(x)  cos  (ut  +  (14) 


-10- 


C2=  -AER^Ej^  CCx  -  En  a  SSx  *  E^t  CCx  +  E2B 


SSx) 


+I^(E1aCCx  +  E^  BSSx  -  E2S  CCx  +  E^  a  SSx) 


+R2(-E^.SCx  +  E  gCSx  +  E2B  SCx  +  E2  aCSx) 


+I2(EjlB  SCx  +  E^  a CSx  +  E2  aCSx  -  E2  BCSx) 


The  free  body  diagram  of  a  package  is  illustrated  in  Fig.  3.  The  equilibrium 
equation  of  "he  package  is 


W°  +  'W\+i>’WIV  Cdpk  >“  *  P"K  at“  =  C 


(17) 


K+l  K+l  K+l 

xK~  o 


m>k  dtlatff  r 


*♦« 


th 


where 


Fig.  3.  Free  Body  of  the  k"  Package 
4owcukakua“ 


'DPk 


as  derived  in  (B.23) 


th 


CQk  is  the  dimensionless  drag  coefficient  of  the  k  package. 


th 


PH,  is  the  mass  of  the  k  package  plus  the  virtual  mass.  The  boundary 

K 


and  continuity  conditions  are: 


(1)  Ui(L1»t )  =  Xq  cosuj  t 


(^)  ?K1Upi  +  CDpi  Upi  -  A1(axl)x]=0  +  A2(ox2)x2=0  =  0 

(3)  U1(0>t)  =  Upi(t)  =  U2  (L2,t) 

(4)  PM2Up2  +  CqP2  Up2  -  A2(cx2)x2=0  +  A3(ox?)x3_0  =  0 


(2n-l)  U^O.t)  =  Up  (t)  =  Un(Ln,t) 

n-1 


(18) 


(2n)  U  (0,t)  =  0 
n 


-ia- 
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From  Un(0,t)  =  0,  R2n  =  I2n  =  0 

By  using  Up  (t)  =  U.  (0,t)  and  substituting  (14),  (15)  into  (18),  the  equation 
k  K 

can  be  presented  in  the  matrix  form  shown  on  page  13,  where 
K12  =  sin  a  cosh 

K21  =  COS  aKLK  3inh  yKLK 
K1±  =  sin  «kLk  sinh 

K22  =  COS  aKLK  COsh  V'K 
pll  =  -ak  eik  aK  +  Ar.  E2K  ^ 
p12  =  AK  Elk  6K  +  AK  E2K  aK 

PK  =  _pm  (JJ2 

13  K 


P14  =  CDP  “ 


P15  =  AK+l(El)K+l(aK+lK22  +SK+1K11  )+AK+l(E2 ^K+l^K+l^l  "^K+l1^?  J 

P16  =  AK+l(El)K+l(aK+lKll  _6K+1K22  )+AK+l*E2*K+l*"aK+lK?2  _BK+1K11  * 

P17  =  AK+1(E1)K+1(BK+1K21  ^  K+1K12  )+AK+l(E2)K+l(aK+lK12  +BK+1KJL2  ) 

P18  =  AK+1 ( E1 ^ K+l (a  K+1K2 1  +0K+1K12  )+AK+l(E2 ^K+lfa K+1K12  'SK+1K21  } 

P21  =  “(AK  E1K  BK+  AK  E2K  “ 

P22  =  “  AK  E1K  aK+  AK  E2K  BK 
P23  =  CBFK“ 

P24  =  -PMK  “2 

P25  3  AK+1(E1)K+1(BK+1K22  K+1K11  )+  AK+1(E2)K+1(“  K+lA  +  6K+1K11  * 

P26  =  AK+1(E1)K+1^K+1K22  +8K+1K11  )+  AK+l(E2^K+lfa  K+lA  ~  BK+1K22  5 

P2?  3  AK+1 ( E1  ^ K+l (  ^  K+1K21  ~BK+1K12  )+AK+l(E2  ^K+l^ BK+1K21  “°X+1KI2  ) 

P28  3  AK+1(E1)K+1(B  K+1K21  ^K+A  )+AK+l(E2)K+lfa  K+1K21  +  BK+1K12  ) 


By  solving  (19),  values  of  ,  1^,  Rj,  I2,  ...  ^n-l’  *2n-l  can  ^  otltained 
and  the  displacement,  strain  and  stress  can  be  computed  from  (14),  (15), 
and  (16). 

A  computer  program  to  obtain  the  equivalent  linear  damping  coefficient 
and  solve  Equation  (19),  written  in  Fortran  IV  and  coded  in  CDC  6400,  is  pre¬ 
sented  in  Appendix  D. 

3.3  Dynamic  Mooring  Line  Tension  Under  Random  Waves 

The  random  vibration  solution  to  a  linear  system  has  been  well  developed. 
The  random  response  of  a  non-linear  system  is  usually  treated  by  either  lin¬ 
earizing  the  non-linear  system  or  generating  a  response  history  by  means  of 
a  simulated  random  loading. 

The  linearization  technique  for  random  vibration  is  much  more  difficult 
than  that  for  the  cyclic  vibration,  and  the  accuracy  can  only  be  verified  by 
experiment.  However,  for  reasons  stated  previously,  a  linearized  system  will 
be  employed  in  this  work. 


3.3.1  System  Linearization  for  Random  Analysis 
Equation  (4)  may  be  rewritten  in  the  form 


<*2u  n  3u  Pr  9ZU  ■  r  n 
U  +  i=-) - cRtr-  t - —  +  E  =  0 

E1  3x2  R  3t  E1  at* 


where 

„  3u  CDT7lDpw  ,  3u  |  3u 

n  3t  ae  Ut1  at 

U  is  a  random  variable,  and  E  is  an  error  vector. 

The  linearization  of  Equation  (20)  may  be  accomplished  by  determining  the 
value  of  C_.  which  would  minimize  tne  square  of  the  error  vector  E  and  then 
deleting  E  from  the  equation,  thus 


9{E2) 

SCR 


{  [CRU 


CDT*1I*D*Pw 


|u|  U]  U) 


0 
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where  {•}  represents  the  expectation. 


Assuming  U(t)  is  a  Gaussian  process  with  zero  mean,  ensures  that  U  is  also 
a  Gaussian  process  with  zero  mean  with  density  function 


p.(U)=  -  -----  e 
U  /2n  Oy 


Substituting  Equation  (22)  into  (21)  and  performing  the  integration,  gives 


„  _  CDT1TDPw  VQ_  V  (23) 

'"R  ~  2AEX  tt  0^2 

Equation  (21)  may  also  be  obtained  by  equating  the  average  power  dissipation 
in  both  systems. 

This  discussion  is  so  far  limited  to  a  definite  point  in  the  rope  to  find 
a  value  of  CD  for  the  whole  rope;  an  averaging  process  must  be  carried  out  by 
integrating  over  the  whole  length: 

f„(|U|.U2idX 
R  "  2AE1  r  L  . 

I  (U2}dX 

*  c°t""d-“w  .  A  1^  <24) 

2AE.  'it  rZ> 

b2dx 

Repeating  this  procedure  allows  the  equivalent  linear  damping  coefficient 
for  the  instrument  package  in  (?.’/)  to  be  represented  by 

o  _  Wpw  /8  . 

CDP - 2  7"  -J  (25/ 


SSSfl&3?TO5*^3&* 


3.3.2  Solution  Technique 

The  spectral  approach  is  a  convenient  tool  for  obtaining  the  statistics  of 
the  response  process  of  a  linear  system  under  random  loading.  The  spectral 
approach  employs  the  equation, 

S  («)  =  |H(io)p*  S  (w)  (26) 

y  x 


wnere 


Sy(w)  is  the  power  spectrum  of  the  response  process. 

S  (w)  is  the  power  spectrum  of  the  input  loading  process. 

H  (w)  is  the  frequency  response  function. 

The  computer  programs  to  obtain  the  equivalent  linear  damping  coefficient 

and  the  spectra  of  the  response  quantities,  displacement,  strain  and  stress 

along  the  mooring  line  are  presented  in  Appendix  E. 

The  assumption  that  the  sea  state  is  a  stationary,  Gaussian  process  with 

zero  mean  ensures  that  the  response  processes  are  also  stationary  and  Gaussian 

with  zero  mean.  The  variances  of  the  response  processes  y,i^  and  are^10^ 
*  dt2 

fOO 


‘SO 

-  S  (m)  doo 
.  o  ” 


J.2  =  J  o?S  ;  ojdw 

y  jo  y 

=  [  ■  (u)doj 

y  jo  y 


and  the  average  frequency  of  the  response  process,  f  ,  is'  ' 

G 

i  a- 

f  =  -kJL 

e  2  it  a 

y 

The  distribution  of  the  peak  values,  yp,  has  been  shown^11^  to  be 
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?  (1— G  )2  2 

'  tn  - : -  x* 


)  (n)  =  -i—  ee  ^c  +(l-e2)~ 
n  /ov  - 


a.  _  n  rn  - - 


e  dX 


where 


n  =  y 


?  1  u 
o  zor.  oB 


Oy2Oy?- 


yp  is  the  peak  value  of  y 

t  is  a  tendwidth  indicator.  The  density  curve  pp(n)  is  shown  in  Fig.  4  for 
several  values  of  c.  j>($) 


Fig.  4  Graphs  of  p^(n),  the  probability  distribution  of  the  heights  of 

maxima  (n=  yp/o  )  for  different  values  of  the  width  e  of  the  energy 


spectrum . 


ft 
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CHAPTER  4.  DESIGN  CONSIDERATIONS 


The  results  of  the  previous  work  suggest  the  following  critical  design 


feature;;. 


4.1  Damping  Parameters 

The  damping  sources  of  a  taut  mooring  line  consist  of  the  internal  damping, 
the  tangential  drag  on  the  rope  surface,  and  the  drag  on  the  instrument  pac) age. 
The  relative  significance  of  the  damping  sources  are  shown  in  Figs.7  and  8. 

The  tangential  drag  on  the  rope  surface  is  clearly  a  major  damping  source.  The 
frequency  response  function  is  heavily  dependent  on  the  magnitude  of  the  tan¬ 
gential  drag  coefficient  as  shown  in  Fig.  9.  Evidently,  a  reliable  estimate  of 
the  tangential  drag  coefficient  is  essential  for  an  accurate  assessment  of  the 
dynamic  force  in  the  mooring  line.  Such  reliable  data  is  not  presently  avail¬ 
able  . 

The  selection  of  a  mooring  rope  should  reflect  the  realization  that  differ¬ 
ent  values  of  tangential  drag  coefficients  will  result  in  different  dynamic 
behavior.  The  use  of  a  rope  with  a  very  smooth  surface  may  introduce  a  high 
resonant  peak. 

4 . 2  Dynamic  Force  Attenuation  Along  the  Mooring  Line. 

The  changes  of  the  variance  of  the  dynamic  force  along  the  mooring  line 
are  shown  in  Fig.  10,  and  the  changes  in  the  frequency  response  functions  are 
displayed  in  Figs.  11,  12,  and  13.  It  is  evident  that  the  higher  the  damping, 
the  greater  the  attenuation  of  the  dynamic  force  along  the  mooring  line.  For 
moderate  or  high  damping,  the  maximum  dynamic  force  is  ct  the  buoy;  for  low 
damping,  the  maximum  force  can  be  at  depth. 
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Figure  12.  The  effect  of  frequency  on  the  dynaaic  force  along  the  mooring  line 


-26- 


<t.3  Anchor  Liftini 


The  mooring  line  force  at  the  anchor  results  from  three  sources:  configur¬ 
ation,  wind  and  current  drag,  and  wave  excitation.  The  chance  of  anchor  lifting 
can  he  minimized  by  adjusting  the  nylon  scope  and  the  rope  combination  (if  a 
compound  line  is  used)  so  that  the  dynamic  force  at  the  anchor  is  a  minimum. 
Another  alternative  is  to  insert  special  dampers  on  the  mooring  line  to  cause 
attenuation  of  the  dynamic  motion  for  some  distance  from  the  anchor.  The 
damper  should  be  designed  so  that  the  drag  coefficient  is  large  in  the  axial 
direction  and  small  in  the  lateral  direction. 
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CHAPTER  5.  SUMMARY  AND  CONCLUSIONS 

Computer  programs  for  the  steady  state  mooring  line  tension  due  to  geom¬ 
etry,  gravity,  average  current  and  average  wind,  and  for  the  dynamic  mooring 
line  tension  under  the  action  of  si:.e  wave  or  random  waves  are  presented  in 
this  report.  The  dynamic  program  is  a  frequency  domain  solution  based  on  a 
linearized  structural  system  and  is  only  applicable  to  a  taut  line  mooring. 
Specific  points  in  this  report  are  now  summarized. 

Analytical  Features 

The  non-linearity  of  a  mooring  line  originates  from  three  sources: 
hydrodynamic  drag,  line  curvatures,  and  material  properties.  Solutions  obtained 
in  the  time  domain  by  using  a  simulated  wave  to  generate  the  output  history 
from  either  an  analog  or  a  digital  computer  are  most  satisfactory  for  mooring 
line  dynamics.  The  high  cost  and  special  'Knowledge  required  in  this  method  may 
not  make  it  readily  available  for  general  design  use.  Alternatively,  the 
dynamic  response  of  a  buoy  system  under  random  input  can  be  obtained  with  a 
small  amount  of  computer  time  from  solutions  in  the  frequency  domain  based 
on  a  linearized  structural  system.  This  method  is  developed  here.  The  linear¬ 
ized  structural  system  has  the  following  features: 

(1)  The  line  curvature  non-linearity  is  neglected  by  treating  the 
mooring  line  as  a  straight  string. 

(2)  The  mooring  rope  is  assumed  to  be  made  of  step-wise  linear,  viscous 
material,  and  the  loss  modulus  is  taken  to  be  always  constant. 

(3)  The  hydrodynamic  drag  is  linearized  through  the  principle  of 
equivalent  linearization. 

For  a  taut  mooring  line,  the  first  and  the  second  assumptions  introduce  neg¬ 
ligible  error;  the  third  may  bias  the  solution  significantly  in  a  high  sea 
state.  No  error  bound  is  available,  and  the  solution  can  only  be  verified 
from  exnerimental  data. 
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forcing  Conditions: 


ia 


The  buoy  is  assumed  to  follow  the  sea  surface,  and  the  sea  state  is 
assumed  to  be  a  stationary,  ergodic,  Gaussian  process  with  zero  mean.  Either 
a  fully  developed  sea  or  an  average  sea  spectrum  is  used  as  the  loading  spectrum 
in  the  analysis  input. 

Based  on  the  linearized  structural  system,  the  dynamic  force  can  be  con¬ 
sidered  as  a  stationary  Gaussian  random  process  with  zero  mean,  and  the  vari¬ 
ance  of  the  process  is  then  derived.  The  total  mooring  line  force  is  obtained 
from  the  superposition  of  the  steady  state  force  due  to  current  and  wind  and 
the  dynamic  force  due  to  waves;  the  total  force  is  represented  by  a  density 
curve.  The  distribution  of  the  peak  transient  farce  is  shown  to  approximate 
the  Rayleigh  distribution. 


The  damping  force  on  a  mooring  line  comes  from  three  sources:  internal 
damping,  water  drag  on  the  rope  surface,  and  water  drag  on  the  instrument 
package.  For  a  typical  mooring  without  sub-surface  buoy  or  special  dampers 
inserted  on  the  line, the  water  drag  on  the  rope  surface  has  been  shown  to  be 
predominant.  The  dynamic  behavior  of  the  mooring  line  is  heavily  dependent  on 
the  value  of  the  tangential  drag  coefficient.  Neither  theoretical  solutions 
nor  experimental  data  are  available  for  the  estimation  of  the  tangential  drag 
coefficient  of  a  rope  under  oscillating  motion.  This  must  be  remedied  if  a 
better  assessment  of  the  dynamic  behavior  of  a  mooring  line  is  to  be  obtained. 
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APPENDIX  A. 


HYDRODYNAMIC  DRAG  ON  ROPE  AND  BUOY 


A. 1  Hydrodynamic  Drag  on  the  Mooring  Rope 
A. 1.1  Steady  State  Flow 

The  drag  force  per  unit  length,  F^,  exerted  by  a  fluid  of  mass  density  p, 
flowing  with  uniform  velocity  V  in  a  normally  transverse  direction  to  an 
immersed  circular  cylinder  of  diameter  D,  can  be  expressed  as 


F  =  -  -  C 
DN  2  DN 


P  DV‘ 


(A.l) 


The  dimensionless  parameter 


'DN 


is  largely  independent  of  Reynold's 


Number,  R,  in  the  range  100  $  R  5  5  x  105  as  shown  in  Fig.  A.l.  A  mooring 

line  has  2  x  103<  R  <  2  x  10^13^  for  which  CDN=  1.2  for  a  long,  smooth  cylinder. 

Surface  roughness  and  strumming  may  raise  this  figure  to  1.8.  In  mooring  line 

(2) 

design,  a  value  between  1.2  and  1.8  is  used;  for  this  work,  is  taken  as 
1.5.  In  the  case  where  the  fluid  approaches  the  mooring  line  with  an  inci- 
dence  angle,  the  normal  component  of  the  drag  force  may  be  calculated  by'" J 


Frn  =  I V D  v„2  =rV’  v2sin2» 

FDT  =  7  CDT(”D  VT2  =7CDT°*D  v2cos2* 


(A. 2) 
(A. 3) 


where<}>  is  the  angle  between  flow  direction  and  the  mooring  line. 


Little  data  is  available  concerning  the  longitudinal  drag  on  the  mooring 

.  .  ( 12 ) 
line.  Based  on  the  towing  tests  on  finite  length  stranded  cables,  Podes 

suggested  the  tentative  relation. 


CDT  =  °*02  CDN 


It  was  made  clear  that  "The  coefficient  of  the  tangential  force  has  not 


been  measured  accurately,  but  the  results  with  respect  to  the  tangential 


•? i  '■•  K&*'  -  -■*  * 


w 


coefficient  agree  at  least  in  a  qualitative  way  with  the  result  of  other 
(12) 

experiments."  This  point  is  further  emphasized  when  it  is  noted  that 

incidental  flows,  instead  of  parallel, were  used  in  most  of  the  cases  where  the 
tangential  drag  coefficient  was  measured. 

The  lower  curves  of  rig.A.l  show  the  theoretical  results  of  Reid's 
analysis''  J  for  together  with  some  experimental  data  from  towing  tests  on 
stranded  cables.  The  roughness  parameter  X  is  defined  as  ratio  of  the  equiva¬ 
lent  sand-grain  diameter  of  the  surface  roughness  to  the  radius  of  the  cylinder. 

A . 1 . 2  Tangential  Drag  Coefficient  Under  Sinusoidal  Motion 
The  tangential  hydrodynamic  drag,  which  is  considered  to  be  negligible 
for  the  steady  state  mooring  line  analysis,  may  be  a  major  factor  in  the 
response  of  a  deep  sea  mooring  line  subject  to  oscillatory  longitudinal  motion 
at  an  end.  This  response  is  critical  in  the  design  of  the  structural  system  and 
therefore  the  tangential  drag  must  be  studied  carefully. 

The  laminar  boundary  flow  around  a  smooth  circular  cylinder  under  longi- 

(14) 

tudinal  sinusoidal  motion  may  be  obtained  from  the  Navier-Stokes '  equation 

ill  +  i  31  =  p.  .  bV_.  >  d  (  . 

3r2  r  3r  v  3t  ’  r  ^  2 

and  the  time  dependent  boundary  conditions 


3 
1 

4 

I 

% 

31 


1 


•js 

iS 

1 


V  =  V  cosw  t 
o 


at  r  =  d/2  for  t  >  0 


and  V=0  at  t  =  0  (A. 5) 

(A. 4)  and  (A. 5)  correspond  to  the  heat  conduction  equation  in  a  circular 

cylinder  with  sinusoidal  boundary  conditions.  The  solution  of  the  above  equa- 

...  (15) 

tions  are 


—  =  A  cos  ut  B  sin  wt+  C 
o 

where  A,  B  and  C  are  functionals  of  Bessel  functions.  The  detailed  representation 
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(15) 


The  tangential  drag 


of  A,  B  and  C  can  be  found  in  Car.' slaw  and  Jaeger. 

force  per  unit  length  then  may  be  evaluated  from 

„  ,3V, 

F  =  udi  -  irdpv  (-T-) 

DT  3r  r=d/2 

A  ..  r3A  ,  3B  .  3C-, 

=  »dpvV  [r —  cosu)  t  +  ^—  sin  wt  +  — J 

o  3r  3r  3r  ,  _ , 

r=d/2  (A. 6) 

where  v  is  the  dynamic  viscous  coefficient.  Rather  than  complete  this  numerical 

evaluation,  an  order  of  magnitude  is  obtained  by  determining  the  drag  force  on 

an  infinite  plate  under  sinusoidal  motion.  The  solution  of  the  laminar  boundary 

(14) 

flow  of  an  infinite  plate  under  sinusoidal  motion  is 


V(x.t)  =  V  e 
o 


2v 


<  — 
cos(wt  - J —  x) 


(A. 7) 


where  x  is  the  distance  perpendicular  to  the  plate 


Then 


V^is  the  velocity  amplitude. 


(tt^O  =  V  A—  sm  wt  -  cos  wt; 

3x  .  o  V2v 
x=0 

=  V  sin  ( wt  -  7-  ) 
o  v  4 

Substituting  p  =  1.935^’sec  ,  v=  10.9  x  10  6  ^t,//sec  for  waxer  at 

20°C. ,  the  tangential  drag  force  due  to  skin  friction  is 
Fdt  =  Tdpv  Vq  M  sin  (wt  -  £) 

=  3.30  x  10_3Pti  d  V  AH  sin(wt  -  ■£•) 

o  4 


(A. 8) 


(A. 9) 


From  (A. 6)  and  (A. 9)  it  may  be  concluded  that  for  smooth  cylinder  at 
laminar  flow  condition,  the  tangential  drag  coefficient  is  of  the  order  of 
10  3and  is  not  proportional  to  V2  but  vAA  with  a  phase  difference. 

(14) 

Surface  roughness  results  in  form  drag  rather  than  skin  friction,  and 
the  drag  may  be  considered  as  proportional  to  the  square  of  the  velocity.  For 
the  drag  of  an  oscillatory  tangential  flow  on  a  rough  surface,  no  theory  and 
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experimental  data  are  available. 

It  is  suggested  that  the  theoretical  solution  by  Reid  for  the  case  of 
steady  flow  be  used  as  a  guide  to  assess  the  drag  coefficient.  Values  for  the 
roughness  parameter  >  are  taken  as  1.0,  0.2  and  0.01;  the  associated  tan¬ 
gential  drag  coefficients,  from  Fig.  A.l,  are  0.013,  0.008  and  0.0035  for 
plaited  rope,  braided  rope,  and  Nolaro  respectively.  These  values  are  thought 
to  reflect  the  relative  roughness  of  these  ropes. 


A. 2  Wind  and  Current  Drag  on  the  Surface  Buoy 

A. 2.1  Physical  Description  of  the  Large  Discuss  Buoy 


The  large  discus  buoy  which  was  developed  by  General  Dynamics  will  be  used 

(8) 


for  this  study  on  the  merits  of  its  surface  following  property.  The  buoy 
is  40  feet  in  diamecer  and  7  feet  thick  with  a  flat  deck  and  truncated  cone 
shaped  underside.  The  weight  of  the  buoy,  including  the  ballast,  is  about 
2  x  10s  pounds  and  the  moment  of  inertia  about  the  horizontal  axis  through  the 
center  of  gravity  is  7.5  xlO5  slug-ft2. 

The  dimensions  of  the  buoy  are  shown  in  Figs.  A. 2  and  A. 3. 

4.2.2  Buoy  Wind  Drag 


Based  on  the  results  of  wind  tunnel  tests  on  a  scale  model,  Nath 
suggested  the  wind  drag  force  on  a  40  foot  diameter  discus  buoy  as 


(16,  6) 


W  2 
a 


F  =  140  p  .  =70  p  .  W  2  (lbs) 

w  1  air  2  Kair  a 


(A. 10) 


Introducing 

pair 

Then 

r 
‘  w 

where 

W 

=  0.16  W  2  (lbs) 

3 


(A. 11) 


The  wind  velocity  increases  with  height  above  the  sea.  For  a  conservative 
picture  the  speed  at  the  mast  top  (elevation  44’)  will  be  used  in  this  analysis. 
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Using  the  wind  profile  suggested  by  Davenport v 

W.  Z.  1/ 

—  =  (— )  a 
W2  V 


(A. 12) 


where  a  -  8.5  and  30’  $  Z  i  800'  for  winds  over  the  water, 
then  (W)z=44,  =  (£>.  8'5  (H)z=64,  =  0.957(i/)2=64, 

It  is  to  be  mentioned  that  for  very  stormy  conditions  the  air  contains 
spray,  which  will  increase  the  mass  density  and  decrease  the  wind  speed.  The 
spray  can  be  considered  to  be  the  result  of  energy  transfer  from  wind  to  the 
ocean  surface,  and  the  total  momentum  flux  of  the  wind  layer  near  the  sea 
level  may  be  reduced  due  to  the  energy  dissipation  in  the  spray  generation. 
Therefore,  the  use  of  a  uniform  wind  velocity  based  on  the  magnitude  at  the  top 
of  the  structure  may  result  in  a  safe  design  for  the  wind  force,  including  the 
effect  of  the  spray.  However,  little  is  known  about  the  true  effect  of  the 
spray . 

A. 2. 3  Buoy  Hydrodynamic  Drag 

Hydrodynamic  characteristics  of  various  buoy  hulls  can  be  found  in 

Hoeraer/18^  Paquette  and  Henderson^  and  Mercier. ^19"^  The  horizontal  current 

(6) 

frag  on  the  large  discus  buoy  was  suggested  by  Nath  to  be 

V  2 

fdc  *  °-035  *  £<40)2  -f-  p» 

=  22  p  V  2  =44  V  2  (lbs) 
w  c  c 

where  V  is  the  surface  current  velocity  in  ft/sec. 


APPENDIX  B. 


DISCUSSION  OF  THE  MODEL  SIMPLIFICATIONS 
B.I  Surface  Following  Property 

Under  wave  excitation,  the  forces  acting  on  a  free  floating  buoy  are 
buoyancy,  initial  force.,  and  hydrodynamic  damping.  If  the  increment  of  buoy¬ 
ant  force  is  much  bigger  than  the  increment  of  initial  force  under  water  undula¬ 
tion,  the  buoy  will  follow  the  water  surface  closely,  and  the  surface  follower 
buoy  is  named  as  a  result  of  this  phenomon.  The  huge  discus  buoy  is  a  typical 
example  of  this  type.  On  the  other  extreme,  if  the  increment  of  buoyant  force 
is  much  smaller  than  the  increment  of  initial  force,  then  the  buoy  will  not  be 
disturbed  significantly  by  the  wave  and  will  thus  remain  stationary.  The  spar 
buoy  is  a  typical  example  of  the  stationary  type.  Other  types  of  buoys  fall 
between  these  two  extremes. 

The  wave  spectrum  measured  from  a  free  floating  wave  meter  (FFWM)  when 
compared  to  that  from  a  slack  moored  large  discus  bucy  by  Gaul  and  Brown, 
indicates  that  the  discus  buoy  behaves  much  as  the  FFWM  at  frequencies  up  to 


Comparison  of  these  spectra  is  sho'/n  in  Fig.  B.I.  The  increment  of  moor¬ 
ing  line  tension  may  restrain  the  buoy  motior  and  hence  distort  the  response 
spectrum.  However,  the  increment  of  mooring  line  tension  is  usually  small 
compared  to  the  increment  of  buoyant  force  for  the  discus  buoy. 

Fig.  B.I  supports  the  assumption  that  the  buoy  response  spertrum  can  be 
considered  the  same  as  the  wave  spectrum  for  design  purposes. 

It  is  suggested  that  wave  spectrum  be  truncated  at  a  frequency  of  0.35 
for  the  large  discus  buoy.  For  smaller  surface  follower  buoys,  the  effective 
frequency  range  may  be  extended. 

The  computer  program  developed  for  the  dynamic  analysis  of  the  mooring 
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Figure  B.i  Comparative  Wave  Spectra  from  FFWM  and  Monster  Buoy. 


line  is  based  on  the  assumption  that  the  buoy  is  a  perfect  surface  follower. 

For  buoys  other  than  this  type,  a  transfer  function  between  the  wave  spectrum 
and  the  buoy  response  spectrum  has  to  be  established  before  a  good  result  may 
be  expected.  However,  the  assumption  will  provide  an  upper  bound  solution  to 
the  dynamic  stress  of  the  mooring  line  provided  the  natural  frequency  of  the 
heave  motion  of  the  buoy  is  far  from  the  effective  wave  .  frequency.  In  the 
case  of  a  stationary  buoy,  the  dynamic  force  may  be  considered  to  be  negligible. 
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B.2  Straight  String  and  Horizontal  Buoy  Motion 

The  mooring  line  does  not  remain  straight  under  the  action  of  ocean  current. 
The  sag  will  depend  on  the  magnitude  of  mooring  line  tension  and  the  current 
velocity.  The  effect  of  lateral  loads  on  the  dynamic  force  of  a  vibrating 
stretched  string  is  investigated  here. 

Consider  a  stretched  string  with  normal  lateral  loading  q(x)  as  shown  in 


Fig.  B.2. 
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Figure  B.2  Vibrating  Stretched  String  Under  Lateral  Loading. 


From  equilibrium  consideration  in  the  x  and  y  directions, 

Ao  cos<!>-A(c  +  AS)cos{^+  f^-  AS)+  mASU  +q(x)AS  sin<J>  =  0 

,  X  X  do  ob 

ana 

Aoxsin$-A(ax+  r  AS)sin($+  AS)-  mASV  -q(x)A3  cosij>  =  0 


(B.l) 


(B.2) 


For  a  prestretched  string,  the  longitudinal  displacement  U  can  be  divided 

into  a  static  component,  U  (x)  =  e  x  and  a  dynamic  component,  Un(x,t), 

s  s  u 

U(x,t)  =  Ug(x)  +  UD(x,t) 
where  is  the  prestretched  strain. 


&*a*X&@2GSK BRSttSS* 


For  a  deep  sea,  taut  mooring  line: 
£  »  Ug(£)  >>  UD 


V  =  U  . 
o  o/. 


eg  is  of  the  order  of  0.1. 


d  D  8  V 

and  are  usually  of  the  same  order  and  much  smaller  than  except 
at  resonance  frequencies.  However,  a  resonance  condition  is  not  likely  to  be 
developed  due  to  the  internal  and  external  dampings.  Based  on  these,  the 
following  approximations  are  made: 

AS 

sin  <J>  =  4> 

cos  <J>  -1 

3lJ  2.  3V  ^  3u 

e=  W  +  f"  (^}  =  3x  ’  £  is  the  axial  strain  of  the  string. 

After  neglecting  the  second  order  terms.  Equation  (B.l)  and  (B.2)  may  be 


rewritten  as: 


92up  1  a2up  ^  q(x) 


AE  3x 


(B.3) 


2  3x  3x 


_D  _32V  _  1_  3 
*  3x2  "  a2  3t2 


(B.4) 


where  a2  =  E/p  ;  p  is  the  mass  density. 

Equations  (B.3)  and  (B.4)  appear  analytically  intractable.  In  order  to  see 
the  effect  of  the  term  on  the  solution  of  Up,  an  approximate  solution 

is  presented. 

The  equations  of  longitudinal  and  transverse  small  amplitude  vibrations  of 
a  stretched  string  without  lateral  loading  as  deduced  from  Equations  (B.3)  and 
(B.4)  are 
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|  32U  32U 

i  D  -  1  °  =  0  (B.5) 

|  3x2  a2  3t2 

|  12L  -i-ifl  =  0  (B.6) 

|  3x2  b2  3t2 

|  ->  T 

|  where  b2  =  jjj-  ■,  in  is  the  effective  mass  per  unit  length. 
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|  T  is  xh^  pretension  in  the  string. 

|  The  solution  of  Equation  (B.5)  with  Loundary  conditions: 

1  U_  =  0  at  x  =  0 

1  U_.  =  U  sin  wt  at  x  =  £ 

k  Do 

1  U 

|  is  U_.  +  - sin  —x  sinu>  t  (B.7) 
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|  The  solution  of  Equation  (B.6)  with  boundary  conditions: 

$ 

1 

*5 

1 

U  V=0  at  x  =  0 

? 

| 

’S 

w 

i  V  =  V  sin  u  t  at  x  =  £ 

|  o 

.3 

3 

1 

!  "  v 

|  V  =  - -  sin  ^x  sin  ut  (B.8) 

£  sin-r- 

t  b 

1 

I 

1 

I  Using  (3.7)  ai.d  (B.8),  then 

i 

1 

f  3V  .  b  32V  ™»  (B.9) 

|  3x2  3X  3  3x2  3X 
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g  where  b  _  /  T 

|  a  /  AE 

1 

% 

%  Introducing  (B.9),  Equation  (B.4)  can  be  reduced  to 

I 

|  _  i  32v  = .  (B.10) 

|  3x2  (l+b/a)b2  3t2  AE 

1 

|  which  lias  a  solution 
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V(x,t)  = 


sm 


- —  sin  x  sin  tut 

tog.  o 

B 


\\°& 


dxdx  +  D^x  +  D2 


(B.ll) 


where 


B2 


=  (1  +  — )b2  ,  and  the  constants  Dx  and  D2  can  be  determined  from 


V=o  at  x  =  0,  t  =  0 

V=0  at  x  =  g,  t  =  0 

Using  Equation  (B.ll)  in  (B.3)  gives 


3X2 


_1 

i2  at2 


q(x) 

AE 


V  £ 

Vo  B 


•  _  to£ 
sm 


cos 


^  »t  -J  agi  dx  + 


For  a  special  case  where  q(x)  =  qQ  =  constant 


D1  =  +  2AE 


d2  =  0 


and  Equation  (B.12)  reduces  to 
32Un 


1  32UD  qo 


3X2 


3t2 


AE 


V  “ 

o  B 


tox 


sin 


.  cos  sinu  t  - 
tog  B 


q  x 
-2_  + 

AE  2AE 


which  has  a  solution 

U.  =  (C.  sin  ^  x  +  C2  cos  ^  x)(C3  sinX  t  +  C4  cosX  t) 
D  1  a  ^ 

^  t,  w 

+  o  B 


(“1  _  ~)AE  sin  ^ 
a2  B2  B 


cos  sinu)  t  - 

B  6A2E2 


q  2  q  * 

—2 —  x^  +  ° 


“TT 
4A  E 


(B.12) 


(B.13) 


(B.14) 


The  last  two  terns  in  (B.14)  involve  only  spatial  variables.  They  may  add  to 
the  statio  component  and  be  neglected  in  the  dynamic  analysis.  With  this  in 

mind  and  letting 

%  Vo  I  CB.15) 


Q  =  - 


(“1  _  “1)AE  sin 
B2  a2  B 
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Then 


Up(x,t)  =  (C^  sin  ~  x  +  c2  cos  ^  X^C3  s^n  **  +  C4  cos  X*) 


+  Q  cos  |j-  x  sin  tot 
Introducing  the  boundary  conditions, 


(1)  uD  =  0 


at  t  =  0 


(B.16) 


(2)  UD  =  0 


at  x  =  0 


(3)  Up  =  Up  sin  wt  at  x  =  £ 


Allow  Equation  (4.30)  to  be  expressed  by 


U 

U  (x,t)  = - 2 - - - 

[sin  —  £-Q(cos  —  - 


- -  sin  —  x  -  Q  cos—  x 

w£v  a  a 

COS  ~g  > 


+  Q  cos  —  suiii)  t 

D 


(B.17 ) 


Comparison  between  Equations (B. 7)  and  (B.17)  indicates  that  the  effect  of 
lateral  loading  on  the  dynamic  force  of  the  string  is  negligible  if  Q«  1.  On 
the  examination  of  (B.16),  it  is  seen  that  Q  is  usually  less  than  unity  by 
several  orders  except  when  sin  0.  This  is  the  condition  where  the  trans¬ 

verse  resonance  occurs,  and  the  solution  based  on  small  amplitude  vibration 
does  not  apply. 

In  the  case  of  a  deep  sea  mooring  line,  the  interaction  between  water  and 
rope  will  damp  the  transverse  motion  and  the  effect  of  the  lateral  loading  may 
be  considered  to  be  small  for  all  conditions. 


B.3  Damping  Linearization 


Based  on  the  principle  of  equivalent  linearization,''  '  the  solution  of  a 


non-linear  equation 

H  X  +  K  X  +  uf(X,X)  =  0 

may  be  approximated  by  the  solution  of  an  equivalent  linear  equation 


(B.18) 
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h 


\,tf 


M  X  +  (1<  +  K^)  X  +  AX  =  0 


(B.19) 


The  equivalent  parameters  and  x  are  obtained  by  equating  the  work  per 


cycle  for  the  two  systems.  The  magnitude  of  the  er.^or  depends  on  the  linear¬ 


ized  quantities  (5-)  and  (4-)  and  is  of  the  order  of  magnitude  of  p2. 


general,  if  pf(X,X)  is  small  compared  to  M  X  and  K  X,  the  error  will  be  small. 


The  parameter  A  is  obtained  by  equating  the  active  components  of  energy  in  both 


systems,  by  equating  the  reactive  components.  Assuming 


X  =  a  cos  (ut  +  <j>) 


(B.20) 


X  and  have  been  shown  to  be 


f 

A  =  _H_  f(a  cos  a  ajsin  41)  sin4>  dij> 

ITclU)  I 
;  0 

r2ll 

K  =  f(a  cos  4, ,—  a  ajsin  4,)  cosij)  d.<j> 

±  7Tcl 

'  n 


(B.21) 


(B.22) 


where 


2_  £ 


The  equivalent  linearized  damping  coefficients  of  the  hydrodynamic  drag 


on  the  instrument  package  and  the  mooring  line  are  derived  as  follow;.. 


(A)  Instrument  Package 


The  hydrodynamic  drag  on  the  instrument  package  is 


FDP  =  ?CDN  ApJ*>  *S»«|X|  X=*X 


where 


CDN  is  the  dimensionless  normal  drag  constant 


A  is  the  projected  area  of  the  package  on  the  plane  perpendicular 


to  the  motion. 


From  (B.21)  and  (B.22) 


X  =  3tt  p“  CDN  A  a  “ 


(B.23) 


Kl=  0 


(B.24) 
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The  magnitude  of  error  will  depend  on  the  weight,  volume,  and  geometry  of 
the  package.  It  will  be  smaller  with  a  slender  current  meter  compared  to  a 
glass  ball.  When  the  water  drag  on  the  packages  are  not  the  major  system 
damping  force,  the  linearization  technique  may  be  applied  to  all  packages 
without  excessive  solution  distortion.  However,  the  arguments  only  provide  a 
guide  for  subsequent  employment  of  engineering  judgment  and  the  accuracy  may 
only  be  checked  by  comparison  with  experimental  results. 

(B)  Mooring  Line 

The  hydrodynamic  drag  on  the  mooring  line  surface  per  unit  length  is 
given  by  F^R  =  ~  A  pjx|  X  =  p|x|  X=A  X  where 

is  the  tangential  drag  coefficient 

A  =  tiD  is  the  surface  area  per  unit  length. 

From  (B.21)  and  (B.22) 

4  a 


Kl  = 0 


(B.26) 


Normally  A=A  (x). 

However  to  ensure  an  analytical  solution  to  the  linearized  equation,  A  has  to 
be  assumed  constant  for  the  whole  section  of  the  mooring  line.  This  equivalent 
coefficient  is  obtained  from  the  balance  of  energy  dissipated  in  the  whole 
section  of  the  mooring  line  in  one  cycle;  this  results  in  the  equation 

fLfT  •  •  fLfT  CDT  • 

(AX)  X  dt  dX  =  ttD  |X|  X)  dt  dX 

•  A  *  A 


giving 


w  La3(x,dx 


a2(X)dX 


(B.27) 


rr- 


The  tangential  drag  is  usually  very  small  compared  to  the  inital  force 
and  the  elastic  restoring  force;  because  of  this,  the  method- of  equivalent  lin¬ 
earization  usually  achieves  good  results. 

B.4  Internal  Linear  Damping  of  the  Mooring  Line 

The  areas  of  the  hysteresis  loops  under  cyclic  axial  tension  of  the  nylon 
ropes  have  shown  that  the  internal  damping  is  not  linear.  When  the  internal 
damping  is  a  major  damping  source,  then  the  non-linear  damping  may  be  linear¬ 
ized  by  using  a  stress  dependent  loss  modulus.  The  problem  may  then  be  solved 
by  an  iteration  process.  Usually  the  effect  of  the  internal  damping  is  rel¬ 
atively  small  compared  to  the  hydrodynamic  damping;  therefore,  only  a  rough 
estimate  of  the  loss  modulus  is  used  in  the  solution  program. 

B.5  Strumming 

Little  is  known  about  mooring  line  strumming  and  its  effect.  The  purpoc'. 
of  the  assumption  of  neglecting  the  effect  of  strumming  is  to  simplify  the 
problem  so  that  it  can  be  handled  analytically.  However,  the  effect  of  the 
small  amplitude  transverse,  motion  has  been  shown  to  be  negligible  in  the  mooring 
line  dynamics,  and  it  may  be  expected  that  the  strumming  will  stay  in  the  small 
amplitude  region  because  of  the  damping  forces  and  the  irregular  profile  of 
the  ocean  current. 
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APPENDIX  C.-  Program  STEADY 


C . 1  Purpose 


To  solve  the  steady  state  mooring  line  tension  under  the  action  of  wind 


and  current. 


C .  2  Program  Logic 


The  Program  starts  with  the  estimation  of  the  inclination  angle  of  the  top 
segment.  Knowing  this  and  the  horizontal  drag  coefficients  in  the  buoy  and  the 
mooring  line,  the  force  of  the  top  segment  can  be  calculated  and  the  position 
of  the  end  of  the  segment  is  determined  given  the  stress-strain  relationship 
of  the  mooring  line.  According  to  (1)  and  (2),  the  force  increment  AT  and  the 
angle  increment  A<{>,  then  the  end  position  of  the  second  segment  can  be 
determined.  In  the  same  manner,  the  forces  and  end  positions  of  the  successive 
segments  are  determined.  If  the  vertical  position  of  the  anchor  thus  obtained 
does  not  fall  close  to  the  water  depth  within  a  pre-set  tolerable  limit,  then 
the  inclination  angle  of  the  first  segment  is  revised  and  the  calculation  is 
repeated.  The  iteration  process  repeats  until  a  satisfactory  solution  is 
achieved. 

C.3  Notation 

(A)  Input 

NLK  Number  of  cases 

K  Number  of  line  sections 

CDS  Current  drag  coefficient  on  the  buoy;  the  drag  force 
=  CDB*  Vc2;  Vc  is  current  velocity. 

CDW  Wind  drag  coefficient  on  the  buoy;  the  drag  force 
=  CDW’  Vw2;  Vw  is  wind  velocity 

CDT  Tangential  hydrodynamic  drag  coefficient  on  the  rope;  the 

1  O 

drag  force  =  ^  *CDT  *  it  •  D’V^  per  unit  length. 
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CDN  Normal  hydrodynamic  drag  coefficient  on  the  rope;  the 
drag  force  =  jpw  CDN.  D  .V^2  per  unit  length. 

VCO  Surface  velocity  of  the  ocean  current  with  zero  wind 
velocity,  in  knots. 

LP  Length  of  the  package,  in  feet. 

WP  Weight  of  the  package,  in  pounds. 

CDNP  Normal  hydrodynamic  drag  coefficient  on  the  package;  the 
dra«f  force  =  CDNP  V^2;  lb. -sec2 /ft. 2 

CDTP  Tangential  hydrodynamic  drag  coefficient  on  the  package; 
the  diag  force  =  CDTP  V^2;  lb.-sec2/ft. 2 

SA  Length  in  feet  oF  the  small  increment  of  the  rope  considered 
to  be  straight.  It  must  be  a  common  factor  to  all  line 
sections. 

DEPTH  Water  depth  in  feet 

ERR  Tolerable  error.  The  allowable  difference  between  the  water 
depth  and  the  anchor  depth  is  ERR  X  DEPTH. 

HS  Significant  wave  height  in  feet. 

UFA  Average  ultimate  strength  of  the  synthetic  ropes  in  pounds. 

L  Length  of  the  rope  in  feet 

DIA  Diameter  of  rope  in  inches 

UF  Ultimate  strength  of  the  rope  in  pounds. 

UWT  Weight  of  the  rope  in  water  in  plf. 

ZPERM  Permanent  strain  of  the  rope  at  station  (due  to  the  anchor 
weight) 

HC  Material  code:  MC=0  for  fiber  rope,  MC=1  for  wire  rope 

NN  Interval  of  the  line  increments  for  output  printing. 
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(B)  Output 

X,Y  Horizontal  and  vertical  axis  with  the  origin  at  the  buoy 

in  feet. 

C . 4  Remarks 

(1)  The  stress-strain  relationship  used  in  this  program  is: 

Z  =  0.171  Exp(-  +  ZPERM  for  the  nylon  line 

and  Z  =  0.0115r  +  ZPERM  for  the  wire  rope 

where  Z  is  the  strain 

r  is  the  ratio  of  the  force  to  the  dry  ultimate  strength.  The  program 
should  be  modified  for  other  stress-strain  relationships. 

(2)  The  convergence  rate  of  the  solution  depends  on  the  stress-strain 
relationship  and  the  initial  angle  correction  function.  If  the  stress-strain 
relationship  is  changed,  it  is  desirable  to  revise  the  correction  function  for 
a  faster  rate  of  convergence. 

(3)  The  program  capacity  is  limited  to  20  line  sections  and  1000  incre¬ 
ments.  The  care  storage  needed  for  this  program  is  40  K.  The  number  of  line 
sections  and  increments  can  be  enlarged  by  changing  the  program  dimensions. 

(4)  For  a  typical  deep  sea  mooring  with  300  line  increments,  it  takes 
3  or  4  seconds  in  CDC  6400  to  obtain  the  results. 

(5)  Output  information  is  self-explanatory. 

C.5  Program  Listing 
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1jiwG«AM  jltADY  (  IimHJI  »oUlHuI  *1Af'E3=I,',Juf  ♦TA*Jt:6='-'d7iJU7  > 

D I  ••«£«£>  I  Ui,  r'hi  ( luo2  )  *(  llvW)»A(lv>i/l)»Y(ioUl)»n{loul)»A(lt-ul)  20  ) 

)  »i-  (  2o  )  »UF  (  20  )  »DIA  (  20  J  »UWT  (  20  )  »^p  <  20)  »MP  (21  )  *i»T(  1001)  *!-^<20>  *CD<»P  <  2 
) >, )  ,CDTP(  ?-  )  20  ) 

REAL  L»LL»LP*LG 

10  FORMAT  (  2  •'  1 4 ) 

11  FORMAT  ( 5F  8 • o ♦ 14 ) 

12  FORMAT  (  4F8.0) 

13  FORMAT  14F6.0) 

14  FORMAT  (  5  F‘  6  •  o ) 

15  FORMAT  ( 2 F 8 • o ) 

READ  lv>*  ‘*LK.» K 
PRINT  21*  NLK. * K 

21  FORMAT  {*  nLA  =  *  *  I4»*  k  =  *jI4) 

READ  14,  CDB.CDw  *CDT,CDR*VC'-' 

PRINT  24,  CDB»CDW,CDT  ,CDiM»VCu 

24  FukmAT  (*  CDB=*»Elu.3»*  CDW=*  »E10  •  3  »*  CDT=* »E10. 3 »* 

1  CDN=* » El-1  •  3 » *  VC0=*»E1^«3) 

READ  13,  5 LP( I ) ,rtH( I ) ,CUNR( I ) ,CDTP( 1 ) »  I=1*M 

PRINT  23 

p?  r  K.'.AT  (*  LP(1)  *,*WP(1)  *,*CDNR(I)  *»*CDTP(I)  **,f 

II*) 

PRINT  33*  ( LP ( I ) , WR l 1 ) ,CUhP ( 1 ) ,CDTP( I ) » 1 ♦  I=1»M 
33  FORMAT  (4Flu«5»IlU  ) 

DO  2 ooO  ML=1»NL<. 

KEAO  12*  jA»DEPTH,ERR 
PRINT  22*  OA, DEPTH, LKR 

22  FORMAT  ( 1H1 »  *  GA  =  * , Fo , 2  * “  D£PlH=* *F 7, 1 »*  ERK=*»T6.:>) 

READ  15*  MO, UFA 

PRINT  2f'»  HS.UFA 

25  FORMAT  (*  HG=*»F4. 1 » *  UFA=**F8.2) 

t'EAD  11*  (L«  1  )  ,DI  A(  I  )  ,i'F  (  I  )  »U.<]  (  1  )  »2PErt',,l  I  )  I  )  »  I  =  1»M 
PRINT  28 

28  FORMAT  (*  U(I)  *»*  DlA(I)  *»*  uF  (  1  )  x»*  0*' 1 

1(1)  * ,  *  2PErm(1)  *»*  Mv_(i)  *»*  1*) 

Pi\IrT  38,  (l(  I  )  » D I  A  (  I  )  *ur  l  1  )  t'O.-.T  T  I  )  »2PERM  I  )  *  »'*  C «  1  )  »I  »  I=1»n) 

38  T  ORMAT  (  5E 14. 3  » 2  I  lo  ) 

READ  iC»NN 
DO  45  I =1 »K 
DIAT  I  )  =L>  I A  (  I  )  /  1 2  « 

45  CONTINUE 
MP(  1  )  = 1 
DO  5-  1=1, K 

MP ( I + 1 ) =MP ( I ) +1 F 1 X ( L ( I ) /OA )  +1 
5 0  CON T I N JE 

C  CALCULATE  INITIAL  ANGLl 
LL  =  • 

DO  3~-  J=1»K 

IF  (  MC  ( J ) • EG • 1 )  GO  Tu  3^1/ 

lL  =  LL  +  L(^)“(  1  ,  +  2Pl><:*‘(  )  ) 

3  COM  TRUE 

l  • 

DO  31-  J=1»K 

IF  ( MC ( J ) • E 0 . 0 )  GO  TO  3iu 
L  ~>  =  L  o  +  L ( J) < ! 1.  +  2PERMI J) ) 

31-  CONTINUE 
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mmuz 


.■•a=®C^wsfiR«K5S!i5i^ 


..~LLf*_G 

.,f_uPL  =  LL/  1  l-LP  !  H~L o  ) 

ZA-iDLPM  •>  ’  /LL 
Ll~>.  171  •'•'/. 

ZL  =  Al_OG(ZZ  ) 

K  A  =  -  • '-  8  1  9  /  Z  L 
F  V=KA*UFA 

ONi"  K MU 1  =  1  •  688 87  FP-> 

IF  { l|j •  Gi_ •  2 4 •  )  00  TO  i2*j 

V/<^12.53R-vwi'T  (HU)  *u«.957 
GO  TO  336 

5?6  VW  =  2.'>’fc»(iiO  +  'j.8)’<'0.9b7 
3  36  VC 3=  (  VCO+.v.23*VW)  *1.68889 

HF=CDvvKVa -ABo  l  Va  )  +  CDB*VCb*AijG(VC3) 

R  t 1 ) =^. 

Z (  1  )=o. 

A  (  1  )  —  y/  » 

Y( 1)=-. 
j! ( i ) =  u • 

F (  1)=-. 

N  =  o 

PHI  (  1  )=0. 

<  DD  =  -v  • 

'nl  (  2  )  =AUNUlt  /TV! 
ia  2  r  i  I*  l  w  I  O  P  =  1 

1^./  F  M  —  1  If* 

Pi\  1  i>t  I  111?  Phi  (  2  )  *TH*  VDt)»Y  (n2)  ♦  ioTuH*i'i2 
111  F  w.<i-!A  T  (  A  £  1  5  •  5  » 2  I A  ) 

FV  =  FH/T  AM<PHI  (?)  ) 

F ( 2 ) =G0RT ( FH*FH+FV*FV ) 

KD=- 

DO  1-,-U  J  =  1 »  K 
:-i  1  =  i-'.F  (  J  )  +1 
M2  =  i*iP(  J  +  l  )-l 
93^  U'J  9  33  I--X1»M2 

i<  (  I  )  =  F  (  I  )  /  0  F  ( J  ) 

R0  =  .r''(  I  ) 

IF  I  MCI  J)  .CQ.I  )  Gvj  TO  950 

o  I  K C  j G ~ T  K A  i •'<  •<£ L A T  I  oijGH I  K  I i>YL\ji'i  i\Ok£ 

Z  (  1  )  =  vd  71 « CAP  l*i-.^l';/.0) 

Zu=Z  (  I  )  +ZPEi?K.{J) 

GO  10  955 

G I  !<£bb-->  I  l<A  I  N  i\£LA  T  I  uNbH  1  f  Fu«\  *  Mt  kuPL 

95-.  Z  (  I  )  =v.ull5*K(  I  )  i  Zi3  =  Z  (  I  1  +ZP L^.MJ) 

955  ->  I  iTP=  5 1 W  (  PH  1(1)) 

Coop  =  co j ( ph 1 ( i ) ) 

0 B  =  0 A  * ( I  •  +ZU ) 

>T(I)=3T(I-1I  +  3B 
Gt’.X  =  3(3*  3 1  up 
jtiY  =  3R*C0->P 
X(  i )=X( 1-1 1+06X 
Y  ;  1  ) =  Y (  I  - 1 )+00Y 
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CUKKEN1  PROFILE 

Y  N  =  Y  (  I  )-.5*SB*C0bP 
YM-YN/3.26-1U. 

IF  (YM.GE.O.)  GO  10  513 
VOVCS 
GO  10  51b 
513  VC=VCo/YM**.4 
515  vCN=VC*C0oP 
VCT=VC*bINP 
D=DIA<  J) /GLUT  t  1  ,  +  Z[3) 
FCDN=oR*D*CDN*VCN*ABGl  VCM 
FCD T  =GB*3.1416*D-'CDT*VC1*AB~>1  VCT  ) 

Ft  I  +  1)=F(  I  )+FC1)T-SAsU/.T(J)‘!C'JG^ 

FAV=(T ( I )+F ( 1+1 ) ) /2. 

PHI  (  1+1)  =HHI  (  I  }  +  (  FCD^  +  oA*  UW  I  (  ~» )  *G  1  )  /  F  A7 

935  CONTINUE 
90u  I=i1P(J+l) 

IK  I  )=u. 

Z( I )=^. 

01 NP=3IN(PHI ( I ) i 
CUoP  =  C0 j ( PH  I ( I )  ) 

X  (  I  )  =X  (  1-1  )+LP(  J)*i>Ii\.P 

y  1 1 1 - y i i-i )+lp( J) “Coop 
jT(1)*GT(1-1)+LPCJ) 

Y i\  =  Y  (  I  )  -  •  5*b0*CoSP 
YM=YN/3,2B-10. 

IF  (  YM,GE •  )  GO  TO  613 

V C=VCb 
GO  TO  615 
613  VC  =  VC->/YM**.4 
615  YCN=VC*C0-»P 
VCT=VC*SIKP 

FCDNsCONP ( J ) *VCN* Abb ( VCN ) 

FCDT=CDTPl  J)“v/CT*ABb(VCT  ) 

Ft I+1)=F( i )+FCDT-WP( J)*CobP 
FAV-tFt I )+F ( 1+1) ) / 2 • 

PHI  t  1  +  1 )  =PHI  (  I  )  +  (  FCDN+wP(  J  )  *oI.mkj/FAV 
lulu,  CONTINUE 
N  =  N+i 


INITIAL  mNGLL  CORRECT  I Um 


YD=Y { M2) “DEPTH 

ydd=yd/depth 

Y  0A  =  ABS ( YOD ) 

IF  (YDA.LC.  ERR)  GO  TO  1«1 
IF  (N.GE.io)  GO  TO  1^1 
YE=i.-YDA**G.25 
/DC=YDA**YE 

PHH2)=PHI  (  2  )  *  (  i.+SIGNt  YDC»Y00)  ) 

GO  TO  1J0 

lul  IF  (<D.EQ.o;  ISTOP=X.2 

PRINT  2U»  .\»YO,bCOPE»DEPIri 

FoK.-iAT  t  lit  1 »  *  i  H I  b  Io  THE  i  o  AT  ITEt'ATI^''  <'iU»''BEf<=  x  »  1 4  ?  * 
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iY-0t»MM=**EU.4.h  L1UL  jO^l-’.FO.A,*  MILK  DEHTH  =  *»F».^//// 
19X»  *X*  .  14X»*Y'  .9X  ,*£l^.  Lt'volrt-  »6X  »  *F  uKCr.*  ♦  8X  ♦  *Fu*Ct-/U.  F  •*  >  3X>  Lu 
lu.l0Ari0us»5X*-AiMoLL  lu  Y**9X»*i*) 

r  KI.4T  3o.  (XI  I)  W(  1  1  »jT  (  1  )  »F(  1  )  »K(  1  )  »Z(  1 )  »PriI  $  I)  »  1 »  I=2»  ISTJP.^A) 

K«]HT  3,;,  x  ( iv*2  1  .YJM2)  *ol  ( >'‘2  )  *F(«2)  fKC*'2)  *4(^2)  •»'HI(|-12)  *«2 

3i/  F UKMA  1  l  /£i5.4tll0) 

,  COiMflf-lUt 
,luP 
EHl) 


fxamplf  uf  i)  a  f  a  input 


A  A  * 

u,  16  yJ  •  ~  U  8 

1.3 

H.b 

3u  •  2  •  1 

.264 

4„5 

3  v-  •  •  1 

.264 

u« 

0  •  0  . 

0  . 

3  U  • 

1 0000. 

.  0>,  2 

i. 

530  ,  0 • 

3-'. 

1.5 

5  3  ‘ .:  v,  ^  . 

1850. 

1.5 

5  3  0  0  j  . 

210. 

1.  J 

53  00V/. 

u  .  06  3 
^  •  Oo  3 
u.0'6  3 
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APPENDIX  D.  Program  DYNSIN 


D.l  Purpose 


To  obtain  dynamic  mooring  line  tension  under  sine  wave. 

D.2  Program  Logic 

The  program  starts  with  a  set  of  assumed  equivalent  linear  damping  coeffi¬ 
cients  to  solve  (19)  and  obtain  the  response  quantities.  The  corrected 
damping  coefficients  are  calculated  from  (B.25)  and  (b.27)  and  then  compared 
with  the  assumed  values.  If  the  differences  are  not  within  the  allowable 
limit,  the  new  damping  coefficients  are  introduced  and  the  program  repeated. 

The  convergence  is  ensured  because  a  higher  damping  will  result  in  a  lower 
velocity  response;  and  a  lower  damping,  a  higher  velocity  response.  The 
coefficients  converge  rapidly. 


D.3  Notations 

(A)  Input 

NLK  Number  of  cases 

K  Number  of  line  sectins 

E  Material  elastic  constant  (see  Table  1.)  in  psf. 

EO  Material  elastic  constant  (see  Table  1.)  in  psf. 

Q  Material  visco-elastic  constant  (see  Table  1.)  in  foot- 
pound-second  system 

L  Length  of  the  line  section  in  feet. 

DI  Diameter  of  the  rope  in  inches 
Z  Mooring  line  strain  under  steady  state  tension. 

RO  Mass  of  the  rope  (including  the  entrained  water)  per  unit 
length  in  slug/ft. 

FM  Mass  of  the  package  (including  the  virtual  mass)  in  slug 
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d 


(B) 


LP 

CD 


MC 


MM 

HS 


Length  of  the  package  in  feet 

Hydrodynamic  drag  of  the  package  acting  in  the  axial 
direction 

Material  code:  MC  =  0  for  fiber  rope;  MC  =  1  for  wire 
rope 

Material  model  code:  MM  =  1,  2,  3,  4  (see  Table  1.) 
Significant  wave  height  in  feet 
XINCN  The  length  of  the  nylon  line  in  which  the  variation  of  the 
velocity  amplitude  can  be  approximated  by  linear  rule. 
XINCS  Same  as  above  for  wire  rope 

CDT  Tangential  hydrodynamic  drag  coefficient  on  the  rope 
surface 

UAP  Assumed  linearization  factor  of  the  water  drag  on  the 
package 

DCR  Assumed  liner ization  factor  of  the  water  drag  on  the 
rope  surface 

Output 

UA  Displacement  amplitude  in  feet 

SIGMA  Mooring  line  force  amplitude  in  pounds 

LT  Total  length  of  the  mooring  line  from  the  buoy  in  feet 

I  Index 


D.4  Remarks 


(1)  The  number  of  line  sections  is  limited  to  10.  It  can  be  enlarged 
by  changing  the  program  dimensions.  The  program  requires  a  core  capacity 

of  60  K. 

(2)  The  program  handles  four  types  of  material  models  as  listed  in 
Table  1. 
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(3)  Virtual  packages  may  be  inserted  in  the  mooring  line. 

(4)  It  takes  a  few  seconds  to  obtain  the  results  of  a  3-section  moor¬ 
ing  line  with  a  water  depth  of  10,000  feet  (compilation  time 
excluded). 

D.5  Program  Listing 


rt 

?, 


fs 


% 

! 


41 


■g 

£ 
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g 

1 

I 

I 

vs 

s 

1 

1 

I 

I 

I 


D  T-<i)I *« (  I  U  I  *^U  TPUl  » I APL5= InPUI  »TAPfc.o=:ujIPuT) 

D I  MOD  I  OH  A(  38  ♦  J8  )  *L(  10  )  .Ul  l  lo)  tKu(  iu)  »L(  lu)  »^'VK  10)  »  oP ( 10) » 

1AK  (  lo)  ,bb(  10)  »i>Ct  lo)  »CoI  !<■■)  »CC(  10)  *fcM  10)  »W1>P  (10)  *Z(  10)  *  L-  T  <  6 

10)  ,UA(  60)  trtC(  10  )  *B(  4u)  * 0lGwA(6O)  »«-b  ( 10  )  ,LP  (  10  )  ♦*(  10 )  ,CD  (  10  >  >A«-P(  i.o 
1  )  »HF  T  <  lo)  »*W(  10)  »<-E  <  lu)  »AEA(  W)  *  ACM  10 )  »  A«A  ( 10)  *A*B<  10)  »UA^(  10)  » 
1AU(  lw)  »CD''  (  10)  » CDK  (  K  )  *  UA^C  (  1^  )  tUAPb'  lo)  »DC“<  10)  »DC'<C(  10)  *DOM  10) 
1  *  MM  ( 10  )  »Eo(  10  )  (  lo  ) 

i< HAL  L,Ll  ♦LO»LIi'iCt<»LH,Lo 
RtAD  10*  hL< 

1~  FORMAT  (14) 

DO  6DU0  NL  =  ltNLK. 

READ  1 0 ♦  K 

<D  =  K-1  -a  <DD  =  K-2 

NRB=4*K-2 

NCB=NRB 

READ  12»(h(I)»Fv>U)»^(I)»L(I)»DI(l)*4(I)»Ko(I)  » «-'•-*  (  I  )  *  lP  (  I  )  ,  CD  <  I  ) 

1C  (  I  )  *HM(  I  )  »  I=1»N.) 

12  format  t  3F  i  j.^,7F6.w.»2I4) 

PRINT  32 

3?  Eu«.'IAT(  1H1  , 3X,*E(  I  )  *  ♦7X»,‘E0(  I  )*»6X»*M  I  )*»7X»*D1 <  I  )*  »6X,*Z(  I  )X»7X» 
lsi\w(I)S6X**i-(I)  -  *  7X  * «  Hi-»  (  i  )*  ,6X,*i-P(  I  )».6X*WCD(  1  )*»3X»*<''C(  1  )*»*  '*«•' 
1(1)***  1  *  ) 

print  22*(EU)*C0(I)*^(I)»l)HI)*Z(I)»RO(l),L(I)  ,Pn<  I  )  *  LP  (  I  )  » CD  (  I  )  * 
1 M  C ( I ) *MM( I ) » I » 1=1 »K) 

22  FOR.-. AT  (in  tlv,£ll*4»31  j) 

DO  18u  J  =  1 » K 

DI  (  o  )  =DI  (U)/buM  (  l.+Z(  J)  )/12. 

180  CONT INUC 

READ  16,  HS»T,XINCN»XIKC6,CDT 

16  FORMAT  (5F6.u) 

PRINT  29,  Hb, T ,X I NCN ,X I NCb , C‘  T 

29  FORMAT  <*  Hb  =  x»F4.1,1f  T=*,F4.1,*  x  I  =  *  *  F3.0  »*  XI 

1NCS=*,F5#^»*  CD  f  =  *»F5«3 ) 

W=6.2832/T 

W3= 

READ  18,  (UAP ( I ) *1  =  1*  KO ) 

18  FORMAT  ( 1 -F8 •  -  ) 

READ  17,  ( DCR ( I ) >  I=1»K) 

17  FORMAT  ( i^F6 • ^ ) 

READ  10,  NN 

N I  T  =  ^ 

4u00  IF  (NIT.Gl.6)  GO  TO  401^ 

NIT=NIT+1 
DO  lou  1=1 ,NRB 
DO  15‘J  J  =  1  ,NCB 
150  A ( I , J)=u, 

ICO  CONTINUE 

DO  2ou  J=l»< 

IF  (MM(J).EO.l)  GO  TO  111 
IF  ( MM ( J ) • L 0 • 3 )  GO  TO  113 
IF  (MM( J) .Eo.4)  GO  TO  114 
CO  TO  115 
HI  0(J)=0(J)*W 
GO  TO  115 
113  0£0=0(  J)/Ev.(  J) 

0W3  =  1 • +Wb  R oE0**2 
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£(J)=C(J)+EolJ)*v*b*wLu**2/oAo 

u  (  J  )  =  U(  J  )  *w/U.\S 
GO  TO  115 

114  FvJo  =  E^  (  J  )  **2+J<  D  )  **2 
E(J)-E(J)+F.u(J)/EQS 
w(J)  =  E-/(J)**2*G(J  )  /EuG 
ll‘i  AR=. 7854*01  (  J1*DI  <  J) 

Art  =  •  7854*01 ! J)*D1 U) 

AE!- J)=A«  *EU) 

Aw  (  J  )  =  AR*w  (  J  ) 

jh*  (  J  )  =GuK  f  (E(J)/i<w(J)  ) 

WoP  (  j  )  =  A'/oP  (  j  ) 
rt'GPS  =  WSP(  j)**2 
LG(J)=L(J)*C l.+Z(J)  ) 

G£( J)=C( J)/E( J) 

QEo=QE< J)**2 

CDP( J) =.848*C0<  J)*UAP( J) 

CDi<(  J)  =2.6667*01  (  J)  ***DCM  J)  “CoT/AECO) 
wE'a'=U£  (  J  )  *CDR  (  J  )  *w 
GA=  WSPS-QEW 
GAS=GA*GA 

G3=CDR<  J)*W+GE(  J)  *«SPb 

GBo=GBJ*GB 

GC=0QRT ( GAG+GUb  ) 

GD=2.*( l.+GEG) 

ALP(  J  )  =Swirt  1  (  (GA+uO/GD) 

BEl  (  J)=-bwRT(  J-GA+GO/GD) 

ALPL=ALP( J)*LS( J) 

3ETL=BET(J)*LGIJ) 
G1NH=(EXP(BETL)-EXP(-BETL) )/2. 
COSH=(EXP(bETL)+EXP(-BETL)  )/2. 

J  )  =cs  I  N  l  ALPL)*oINH 
GCt J)=GIN( ALPL)*C0GH 
Co ( J ) =C0G l ALPL ) *o I NH 
CC ( J ) =COS ( ALPL ) *CuSr) 

AEA( J)=AE< J)*ALP( J) 

AE8U)=AEU)*BEI  <  J) 

AUA<  J) =Au( J)*ALP( J ) 

AQ6< J)=AG( J)*BFT ( J) 

20C  CONTINUE 

DO  3ow  I  =  1 » K 

j  =  4  k 1-3  $  JJ  =  J+1 

A  (  J  »  J  )  =oC  t  I  )  *  a(  JJ»oJ)=jCU  ) 

300  CONTINUE 

DO  4Cw  I = 1  * K 

j=4#I-?  1  JJ=J-1 

A ( J  *  J J ) =Co ( I )  o  A( JJ*J)="Co{  I  ) 

400  CONI INUE 

IF  (K.D.LE.U)  GO  TO  1^30 

DO  5ow  I =1 »KD 

J=4* I -3  o  JJ=j+2 

A ( J  » J J ) =CC (  I  )  *  A( J+l*JO+l )=CC(  1  ) 

5o0  CONTINUE 

DO  60-  1=1* KL) 

J  J=4 *  I  j>  J  =  JJ-3 

A ( J  *  J J ) =Go (  I )  oA ( J+l »oJ-l ) =~io (  I  ) 


60u  CON I  INUE 

00  7  u v  1=1*  K.0 

J  =  4  «  I  4  l  i,  JJ-J-Z 

A(J»JJ)=~i.  -i*  A(  J+1»-jJ+1  )=~1« 

7 wo  CONTINUE 

DO  1=1^0 

J  =  4  *  I  - 1  i  Jk  =  JU 
10=1+1 

A  (  J*  J-2  )  =/'Ud  (  1  )  -ALA  (  1  ) 

At  J* J-l ) =  AEU l I ) +  AUA ( 1 ) 

A  (  J  *  J  )  =-k..,(  1  )  x.;o 
At  J*J  +  1  )  =-CD^(  I  )*.V 

At  J* J+2  )  =  (  A  E  A  l  1+1  J-A^fjt  I  +  1  )  )  “CCt  I  +  1 )  +  (  ALB  <  1+1 )+A*A<  1  +  1 )  :*SS  (  1  +  1  ) 
At  J»J+3)=(  At' At  1  +  1  )“A^B(  1  +  1 )  )  -v-ot  1+1  )  ~  (  AEb  (  1  +  1  )+A‘-A(  1  +  1 )  )*CCl  1  +  1  ) 
At  J  +  l *J~2  ) =-A<  J* J-l i 
At  J+1»J-1)=A( J»J-2) 

A ( 0  +  1  *  J  ) =~A ( J  »  J+ 1 ) 

At  J+1*J+1)=A(J»J) 

A  (  j  +  l  *  j  +  2  )  =  (ALB  (  i  v)  +A«A<  1^5  )  'Ll.  t  1  ^  )  +  (  Awot  lv  )“ALM  I  u  )  }  (  lw  ) 

A  (  J+ 1  » J  +  3  )  =  (  ALi>  (  I  O )  +  AuA  (  1 0 )  )'*.ol  w )  —  (  Auii  t  I  w  )  —  A  £  A  (  I  u )  )  *CC  l  I  o  ) 

8  .0  CONTINUE 

If:  (KDD.Ll.O)  GO  7  0  I0  3v/ 
og  i  =  i.<Otj 

J- A  x  1  l  1G=H1 

JA= J-l 

A  l  J.' » J+  3)  -  l  ALOt  W  )  +A^A(  H  )  ku(  I  +  (  A^b  (  1  u )  -ALA  t  I  w )  )  *dC  (  1  ^ ) 

At  Jn.J  +  4)  -(AL:i<  lOJ+A^At  I'-O  )*GC(  i  ^  )  —  (  Av-ti  *  I  o  ) -ALA  <  I  o  )  )  *Cb  t  1 «-  ) 

A  t  J  ♦  J  +  3  J  =  (  ~Al  A  (  ;  v. )  +Av«B  (  I  ^  <  )  x  v.  C  (  w )  - 1  AEb  (  i  o  )  +A>-A  (  1  v )  )+$C(  I--') 

A  i  J*  J  +  4  )  =  (  -At  At  W)+Auj>  (  W)  >  »  bC  (  1^)  +  (  AtBt  lo)  +  A-«A  (  I  w )  )  *  Co  (  I 
9)0  CONTINUE 
lo3>,  DC  1  -,rjJ  1=1* NRo 
D  !  1  )  =  ->  • 
l-'3v.  CON  I  INUE 

dll) =HG/ 2 • 

N  =  NRB 
M=1 

I  r>  I  Z  E  =  3  8 
JG  I  ZF.  =  NCb 

CALL  Ii<V!<(A»  i«»l>*  ‘-I » oL  f  L''*-> »  I  i  L  l  »  Jj  i  L  t. ) 

•<.4=4*<.  t,  K  3  =  k4—  1 
b(K.4)=u«  1  d  l  K3  )='»/, 

C  Pi<INT  5-J*  <B(I)*I  =  1»N4) 

3-  EOK.mAI  (12Eiu»2) 

N=  1 

L  7 ( 1 )=o. 

DO  llwu  I=l*< 

>,2X  =  «;« 
u  3^  • 

J  •„  =  v> 

J  K  -  4  *  i  i  o  o  =  J  >>  -  1 

o  I  J K - 2  *  JU  =  JA-3 

C  1  =  ■(  JI  )  *->C  (  I  )+d  t  Jn  )  *CC  t  1  )+d  t  Jil)  »CGt  i  )-i,(  Jo  )*oo(  I  ) 

C2-t  ( JJ)  "CCt  1  )+i)  t  JO  *GM  I  )  tlj(oil)  >GC(  1  ) <  Ji  )*Ci>(  1  ) 
oA(N,-j<;iti  f  C 1  *  C 1  +C2  *C2  ) 

A^o -f  CAt  I  )  *..j1  i  ) 

ACC  =  A. At  I ) *CC (  1 ) 


AbC=ALA  (  I  )*bC<  I  ) 

ACb=AEA(  I ) *Cbl I ) 

0Gb  =  A£l3<  I  )  *bb(  I  ) 

HbC=ACB( 1 ) " oC (  I  > 

BCC=AEB( I ) »CC (  I  ) 
iiC-j  =  ALB  (  I  )  *Cb  (  I  ) 

C4-BI  JI  )  *(  AoS-BCC  M-dIJnIMiIoU  Au,  )+d  (  Jtl )  *  (  Bbb+ACC  )+ol JJ  )  * (  dLj“AoL  ) 
C  3=L  (  JI  )  *  <  ACC+Bob )  +B  (  Jn)  ( dCb-A^C ) +B  (  0|  j  )*  (  bCC-Abb  ) “B ( o J ) * (  BoC+ACo  ) 
Cb=C4-GE(  I  )  * C 3 
C6=C3+GEl 1 ) *C4 
o  I  G'-'A  (  N )  =  oGl<T  (  Cb*C5+C6*C6  ) 

IF  (MC(  1  l.EG.l)  GO  10  1210 
XI  NC  =  XINC4 
GO  TO  122^ 

121^  XI NC=XIWCo 

122^  X  I  NCb  =  X I  NC*  ( 1  •  +  /■  1 I )  J 

IF  ( L( I ) .LE.XINC  1  GO  Tu  l2uo 
JO-  I  F  I X  (  L  <  I  I/XINC) 

DO  115C  J  =  1  ♦  JD 
N  =  N  +  1 

LT(N)=LT(»m-1)+XINCS 
X  =  Lb ( I )-XlhC5*FLOAT ( J) 

A  L  P  X  =  A  L  P  (  1  )*X 
t3E T  X  =  BE  f  (  I  )  *X 

EXB  =  CXP{  BLTX )  EXB‘'J  =  EXH(-BETX) 

oIImHX=(EXl)-EXBN)/2. 

COGHX=(EXu+EX3N)/2. 
bbX=blN( ALPX)*SINHX 
GCX=SIN( ALPX)*C03HX 

CoX=COS( alhxj^sinhx 
CCX=COS ( ALPX ) *C0SHX 

Cl=o(JI)*oCX  +d(JM*CCA  +u(Jn)*CbX 

C2  =  b(JJ)*CCX  +B  l  JX  )  *ii>X  +B(oH)*oCX  -dIJI)*Coa 

UA(N)=bGKl (C1*C1+C2*C2) 

ASS=ALP( I )  *b5X 
ACC=ALP( I )*CCX 
A5C  =  ALP(  I  )*SCX 
AC3  =  ALP (  I  )  *CoX 
Bbo=0ETT I )*3bX 

bg:=bet ( I )*GCX 
BCC=BtT( I )*CCX 
BC3=bET ( I ) *C3X 

C3  =  B(  JI)*<  ACC+t3S3)+B(0(s)  *  ( BC3~AbC  )  +d  ( OH  )  *  (  BCC-Abb  )  ~B  <  o~»  )  *  (  BoOACo  ) 
C4  =  o  (  JI  )  *(  ASS-8CC  )  +B C  JX. ) ^  ( BbC  +  ACo  )  +b (  OH ) *  (  BSS+ACC  )  +B  (  o J  )  *  (  bCb~AjC  ) 
Cb=C4-GE ( I )*C3 
C6=C3+QE( I ) *C4 

o  I  G>v’iA  ( N ) =3ukT  (C5*Cb  +  C6*C6)*AE<  I  > 

uD=(UA(W)-UA04-l)  )/JA(N-l) 

JDb=UD*UD 

OAo  =  UA(-l“l  )*JA('<~1  ) 

U2X=G2X+JAb*XINCb*( 1  .+UD+UD3/ 3 . ) 

03X=U3X+UAS*UA  ( N”1 )  *XINCS*(  1 .  +  1 . 5xUD  +  UDb+UDS*UD/4  .  ) 
lib 0  CONTINUE 
1201,  N  =  N+1 

NJD=N-JD-1 
LI  <N)=L3( I ) +L T ( NJD ) 
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'JO  U  U  O  O1 


UA  ( N )  =  S  jk  I  (  U  (  Jis )  •>  ij  (  Jk  )  +u  (  J  J  )  -  u  l  J  J  )  ) 

UAPC  (  1  )  =UA  ( N) 

UAFo( I )=(oAH( I )~JAHC( I ) )/JA^< 1 ) 

C3=b(JI  )*ALP(  i  )+o(  Jh)*iiLT(  I) 

C4  =  -R  (  J  I  )  K  BE  T  (  J  )+B(  JM)  *ALI'(  1  ) 

Cb=C4-QE  (  1  )  »C3 
C6=C 3+QE  l  I  )*C4 

S  I GMA  (  N )  =  ->uKT  (  Ci>*CS  +  C6*Co  )  *Al  (  I  ) 

L  lnCK=LT  ( •  J )  -L  1  (Im-1  ) 

UD  =  (  UA  (  M  )  "jA  (  N-  1  )  )  /<JA  (  im-  1  ) 

JDo=UD*UD 

UA0=UA  (  M  —  1  )  *UA  ( i*Ni—  1  ) 

J2X=U^X  +  'JAS*LI.MC‘<-(  (  1 . -t-JD+iJLo/ ->.  } 

0  J  X  -  IJ  3  X  +  UA  S  *  U  A  ( N~1  )  *<-  IRC1'*  (  1 .  + 1 . 3  A  UD+UD3+UDS*UD/4  .  ) 

N  =N  +  1 

L  r  (  i'M  )  =  L  T  (w-D+LRi  i  ) 

DCRC<  I  )  =u3X/U2X 

DCi<-(  I  )  =  (uO'C  (  I  )"DC‘<(  I  )  )/UC'(  1  ) 

11. /o  CONTINUE 

C  I  .Ml  26 »  (OCM  I  )  »OC-<C(  I  )  »DCro(  1  )  ,  1i’,m 

26  FORMAT  (9L13.5) 

00  3»*»,  I  =  1  » K 
DuR  ( I ) rOCKC (  1  ) 

3^6^  CON T I NUE 

00  3-;3v  I  =  1 » KO 
JAP ! 1 ) = ( UmP ( 1 )+oAPC (  1  )  ) /2. 
o«AH  U  (  I  )  =UAP  (  I  )  /(.ARC  <  i  )  - 1  - 
3;3o  CONTINUE 

ALLUaAuLl  ERROR  FuR  OCR  I  j  1'LrCLNT 
00  3  .v  b  U  I  =  1  »  k 

IF  (  AdolOv-i'ol  I  J  )  •  G  T  .  ^  •  1  }  O'o  To  4  o  o  U 
3^8o  CONTINUE 

ALlu.-.AdLt  ERi<Gi\  FOR  GAP  IS  PERCENT 
DO  3^9^  1  =  1, KD 

IF  ( Ad- { j (  i  )  ) »GL  •  ^ • 2 )  Co  lo  4uoo 
3 U 9 o  CONTINUE 
4o  1  -j  NH=  N- 1 

PRINT  8o»  (DCMi),  1=1, A) 

8v  F'-'R-'-AT  (in  » /  /  /  /  3X  »  "DC*'  (  1 )  *  »4X  »  K  OCR  (  2  )  *  » 4X  » *00'  (  3  )  * » 4X  ,  *DCM  4  )  *  *  '*X 
1 » "OCR ( 5 ) *  »4X , "DCR ( 6 ) ^  »4X  j^OCr ( 7 ) *  »4X  »>DCR ( 8 ) *  »4X  *  cDCR ( 9 ) * ,4X  *  *OCR i 

1 1  >)*//(  ]  -  f  i .  2 ) ) 

PRINT  84,  (UAP(I),  I=1,R0) 

P,u  F^R'-‘AT  (  lrl  »////:>*  »  *OAP  (  i  )  *  »4A  »  ,,UA'' (  2  )  “  »4A  »  "oAR  <  3  )  X»4X»  *UAP  (  4  )  *  »  4A 
J  ,  U  A  p  (  5)  *  »  4  X  ,  A  A  p  (  6  )  “  ,4^  »  ”  jA1'  (  7)  *  »4X  ,  *  U  AK  (  8  )  x  »4X  ,>0AK  (9)x»4X*xoAr'( 
llw)  -//(  KFlu.2)  ) 

PRINT  6U*.<*NIT 

6w  F  o>'i-iAT  (In  ,///////// -  I.II^  i  ^  J,ji_  uLiUuTi  .,-lTn  uAi-'Plnu  mT 

1  oE',C  7  =  k  ,  Fo  ,  6 ,  ^  *»U  'i3iir  -F  I  T  E  rA  1  I o.,  =  *  ,  I  2/ / / //9X  » *UAX  *  1 2X  ,xi  i  u  •• 

1A*»  14X  ,*u  * ,  i .  >; ,  k  ; «  ) 

|J|'  I N  T  7  9  ,  (  UA  (  I  )  ,  S I  G-*iA  (I),LT(I)»I,  I  =  1  ,nH  ,  mi ) 

/w  FORMAT ( 3E1S.4, I  1^ ) 

6*jou  CONTINUE 
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-.  I  UP 
l' NO 

';l|;"-UI  "“•  . . . . .  i 

LvrUl  VALEuCr  (  i  /  f 

C  "  A)  »  (  IC^Lyn.jCuLU.-i  )  ,  (  A,.,AX  t  |  ,6l vAH) 

1'-  DE.  I  Li<  -1=1,- 
1  OO  t*  -  J  =  i  *  ti 
2-  IP l  VO  I  (  j )  =  v/ 

3v  uu  ]=it/\ 

C 

C  t At<CH  F0>\  HIVOT  LLf.;-iL'NT 

4-  AKAX=^.- 

4i>  DU  I  ■-  -j  J=x,;g 

‘J,/  1F  '  I^IVuf  i  J,-i  ,  6u,lu0.o- 
6-;  bo  U-  6=1  »;j 

!F  1  IH1VuT  <M-1  )  8b,loo.  740 

DO  I  C0LUM=,< 

05  A.'iA/.  =  A  ( J,..  ) 

1  o ./  C_  oN  T  I  UUh 
l'-b  CON!  Ii\U£ 

r  U'  I  ,JI  V'-'T  ( 1  C^L. jr,  )=l^iVv-<T  (  i  U'-'Hi.'i)  +  j. 

C  I  I  I  L,U_HAi'!u  i>Oa  j  1 1/  ruj  Vl\!  i  i  r  ii-  i 

C  1Vv"  L«.C.hLi»1  ui,  DlAb^iWU- 

rirr!.I;W^“1C°LlJh)  14^’  l*o 
i4>.  bcra^xa-DCTt;?;-: 

1  b  -  oo  2  >  -  L  =  i  * rg 
16-*  -  "  A  P  a  A  (  J  N  o  *\  f  l  ) 

17^  A  (  1  !<«->*»  L  )  =A  (]  i<ua  »L  ) 

2  '^  A(  ICULUMn- j  =u.vAH 

I  F  ( 26.8 »  26  - »  2 1  o 

£i  1  DO  ^  v  —  L  =  i  » i-i 

.8/1  Al’sl)  (  J  I.O.,  ) 

2  3’J  til  ]AU,v,L)=(3I  KoL*j;.;,L) 

2  b  ^  0  (  i  C  0  L—  "A  »  L  )  =  o ..  A  P 

26.8  INDEX!  I  ,  J  )  -  J  i(0\-: 

2/-  1  .mDEX(  1  ,4.-)  =jculu;-« 

31  '  ^.I/oT  <  n*«UCwLU*-.#ICwLa-iJ 

3  2-8  Ol  TE  W;-l  =  Dt  f  L  k;-;«P  Wut  (  j  j 

DIVIDE  I'WbT  ,<o>.  by  PJVo, 

3  3-  A  (  I  CwLUi-i  *  I  CoLL.', )  =  1  ,  j 
34-  bo  j  -j ..  L-i  *N 

Yl !  ’ w  >  =A  (  K^LU.- , L  )  /»■  I  -V8 ,  (  J  , 

3bi>  ir  {.•1)38  8,36..  ,36- 
36-  DO  37-  b=i 

3  7-  R  (  ICoLJ.-i,  -  )  =b  (  IC^LU‘y‘tL)/t'j\y^r  /  j  j 

PiSuouL  ww.-l'iyoi  .w..;6 


"  >c»  I- ?»»*-<* 


APPENDIX  E.  Program  DYNRAN 


K.l  Purpose 


To  obtain  the  dynamic  mooring  mooring  line  tension  under  random  waves. 


rogram  Logic 


Same  as  DYNSIN  except  the  corrected  damping  coefficients  are  calculated 
from  (24)  and  (25). 


E.3  Notations 


(A)  Input 


Same  as  in  Appendix  D,  plus 
F  Frequency,  in  cycles/sec. 

KW  Number  of  frequencies  considered 

NSPEC  Wave  spectrum  code: 

NSPEC  =  1  for  Pierson  and  Moskowitz's  spectrum 


NSPEC  -  2  for  Scott's  spectrum 


(B)  Output 


UA  Displacement  amplitude  under  unit  Sine  wave  excitation 

USPEC  (UA)2  *Sjiji(f )  *Af ;  S^(f)  is  the  buoy  motion  spectrum 

EBSILO  Strain  amplitude  under  unit  Sine  wave  excitation 

ESPEC  (EBSILO)2  S..  (f)  Af 
nn 

SIGMA  Force  amplitude  under  unit  Sine  wave  excitation 

TSPEC  (SIGMA)2  S,,  (f)  Af 
hh 

LT  Mooring  line  position  from  the  buoy 

UVAF  Variance  of  the  displacement  c^2;  in  (lbs.)2 

EVAR  Variance  of  the  strain  a  2 

e 

TVAR  Variance  of  the  dynamic  force  q^.2;  in  (lbs.)2 
■« 

TM2  u?  S  (u)d;  Stt(u))  is  the  force  spectrum 

n 


w 

* 

TM4  to4  Stt(u))d  « 

0  L 

TEB  The  band  width  indicator 

I  Position  index 

VSAR  f  Oy2  dX;  V  is  the  velocity 
' o 

VTAR  [  ov3  dX. 

•*  o 

DCRC  Corrected  linearization  factor  of  the  mooring  rope 

UAPC  Corrected  linearization  factor  of  the  package. 

DCRD  DCRD  =  DCRC/DCR  -  1 

UAPD  UAPD  =  UAPC/UAP  -  1 

K 

DCRS  DCRS  =  l  (DCRD)2 


UAPS  UAPS  =  l  (UAPD)2 
1 

E.4  Remarks 

(1),  (2),  (3),  and  (4)  same  as  in  Appendix  D. 

(5)  If  the  linearization  factors  do  not  converge  to  the  acceptable 
values  in  six  iterations,  the  operation  passes  to  another  case. 

(6)  The  frequency  may  be  arbitrarily  spaced. 

E.  5  Program  Listing 


U 


Ks-S  ^‘:v.'-  vra^sw^i;  f 


K«<-><ji\Ai*.  u  <  iiixAii  (  1 1'lr'u  T  ♦  ^0  T^tj  1  *  i  uHl  5  =  l  *”RU  1  * 1  ARc.6=,RUTRUT  ) 

DI.il'a'jlo.M  A(  38*38  )  *t  (  10)  »DI  l  lo)  ,«J|  10)  »L(  10)  *Mim{  10)  »  b^ ( 10 ) * 

1  Af  <  1  •.  )  *Cb<  T  0  )  ♦  br<  l  -  )  *rb  <  lo  )  »CC<  U1 )  ♦  XX  (  1  v  )  *ivbR<  lu  )  »Z  <  10  )  »F  (  20  1  ><-  T  (  6 
1*.  )  *IJA  (  6> )  »'-‘C  (  1  0  )  »Ri  i|.j)  ♦  :>  1C.*'>A  loo)  *i-o(  10)  *  Lp  (  10)  »u(  10)  »CD  (  10)  *A0H(  10 
1 )  »  Rr  T  (  1  -  )  *v.-i  Uw)  *OF  (  lw)  »  ATA<  1 J)  *  A  FI  B  (  lo  )  »AUA(  lo  )  >Av[}(  10)  *CA F  <  1 J)  » 

)A  1  -.)  »rr>r(  ,r;iK(  lv)  ,UARCt  1-,  *  DC'^t  10)  » Eo  i  10)  101  *V-R 

1EC  (  6*< )  ,VVV'  ( 6„  )  »  ,Ab(  2-j)  »OAH0  (  1  u  )  *L>C‘'C(  lo )  *DCKB(  10  )  »VbA'M  10)  » V  I  AM  1 
).})  * >-E  (  To)  *M  (  lo)  *  TbM  r(  6u)  ♦ob^fiC  (  60  )  *  F.  [io  I  L-<  {  60  )  »EbpEC(60)  *UV  AMoO) 
)  *FVM<(  6  / )  *  TVA*  (  60  i  * 1  •■>2  (  6j  )  ,  f  •-•/*  1 6u  )  *  T  EB ( 60 ) 
a  L  Al  L  * L  I  » !  U  *  L  I  ■  *CR  ♦  lp  *  L  b 
1  M  FOR:' A  f  U  M  l  2.4  *  I  lb  ) 

RFAD  I  '• » NLK  *7. 

1  ,  FORMAT  ( ? I  4 ) 

:yj  c,  i  o  <l  =  i»nlx 

i\  FA  l.>  12*  (:  <  I  )  »Tol  1)  **(  !)  *H  1  >  »l>l  <  1  )  *4(  I  )  » (  I  )  »•'.*.{  I  )  »!-.-(  I  )  »Ci><  i  j  *■• 
1C  (  I  )  I  )  *  I  =  1  »  K  1 
12  FORMAT  i3Fl;.-»/F6*<j*214) 

PRINT  32 

F^.x.-.AK  Dll  *  7X  *  kE (  1  )  -  »7X*KE0l  1  ) x  »  6X**M  I  )*  *7X  **01  l  I  )  *  » 6X  **Z (  I  ) x  *  7X  » 
1  KRu(  !  )  *  »f,X  »*L  l  J  ) *  »7X»*PM  I  l*  ♦  6X»*‘-.-M  I  )*  »6X  >VCD(  1  )*  »3X**,-‘C(  1  )  *  *w  • 
1  i  I  1  *  *  s  I  *  ) 


1  v  1  j  -  7  •'  *  # 

?2  ♦  <  r  n  )  »lo(  i )  *-m  »im  (i )  *Z(  i )  i )  ,n  i )  n  »lh{ 1 ) *cd( i ) » 

1XC  (  1  )  *  M.'-t  l  1  )  » I  ♦  I  =1  »N) 


22  FOR  i  1  1 «  (  » 1  -El  1.4*31  b) 

Dvj  1  hJ  J  =  1*X 

01  (  0)  =DI  (~> )  /bv.l<l  (  1  .  +4  1  J )  )/12. 

1=!.  corn  I 


<0=2-1  1  KDD=l'-2 

RRH=4 ■<-? 

RCB-fiRB 
A' FAD  6  »  X * 

5  FORMAT  (14) 

><  EAD  14*  (Ft  1  )♦!  =  !*  n.-.  ) 

14  FORMAT  (2'/F4.-) 

READ  16*  •  lb » X I  Nf**  *X  3  "*Cb  *CDT  ♦ 

14  FORMAT  (  4^6.'-*  14) 

t'K  [uT  26*  Hb  ♦  X I  MCii » X  I  .<r;>  *C0  l  ♦■'•SPEC 
26  F x/Ri-. A  T  (  s  :  lo  =  *  *  F  0 . 2  *  *  X  1  ««0<s *  »  F  6.  1  *  * 

}  CDT=»*F6.3*S  NSHEC -  * ’ 1 2 ) 

REAu  18*  ( OAR ( I ) *  I  =  1 *KD ) 

18  FORMAT  (  1  -F  6  ••  ) 

READ  17*  <nCHl  1  ) *  1  = 1  *  X ) 

)  -r  FORMAT  (  1  vF6.  • ) 

PRINT  ?7*.s 


XIi’.Cb  =  *  »Fo.  1  *  * 


27  FORMAT  (* 


PRINT  7  7*  OCMD*  1=1*^) 


OCR ( i )  from  i  to  » 


2  8  FOKiAAT  (  1  uAp(I)  Fin~>>'*  l  Tp 

R  a  1  OT  77*  l  u  A  <■'  (  I  )  *  I  =  1  *  nD  ) 

37  FORMAT  ( i - E 1 2 • 3 1 
30P7  1=1*6- 

’JVAR  i  I  )  =b. 
p VAR ( I ) =  b . 

;  VAR  (  1  )  =->. 

T  -’2  (  I  )  =  -  . 

T  M  4  (  I  )  =  —  • 


mi 
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193  COM  I  1  NUL 

IF  ( NSPEC • E w. 2 )  GO  TO  94 
IF  tHG.GF.24.)  H5  =  J  . .  2  7«  (  H8+i> .  a  )  **2 
94  DEC  -  J. 

K  ti  tt  -  K  n  -  1 

DO  128  I  n  =  1  *  K'a  W 
i)F A  =  F  (  IK)  +  F(  1K+1  ) 

DF=DT A/2.-DFC 
DFC=DFC+DF 
W=F ( IK) *6. 2832 

iv  5j  =  a  *  •/. 

DW=DF*6.2832 

IF  tNGPEC.FO.2)  GO  TO  126 

U  A  T  =  (  .  OJb  3d  85  /F  ( I  K )  *  *  6  }  *EXP  (  ~.u2  1  32/H0**2/F  (  In) *#4) 
'JAM  I  ;<  )  =UA  I  *0F 
GO  TO  128 

126  •'<  ’-I  •  /  (  -j  • ^3*H6+ 1  •  35  ) 

W  W  =  •%'  *“A,(/  +  »y  •  2  6> 

I  F  ( '.<V\  L£  •  )  GO  10  127 

WF  =  TJRT  (  (  •  26 ) **  2/t  •065/aW  ) 

U AT  =  t,.214*HG*H$*EXP  (  ~WE  ) 

JAG  1  IK ) =JA1 *DW 
GO  10  128 

127  UAo( IK)=u. 

128  CONTINUE 

PKI.'rr  1  34.  (Ft  j)  .UAo<  J)  ,J=1  ,N/..v  ) 

134  FORMAT  (2F14.4) 

N I  T  =  f 
NKK  =  '» 

2w.yl  00  93  1  =  1  »6o 

WAR  (  I  )  =  0. 

93  CONTINUE 

DO  2  u 0  I  K  =  1  ♦  < >. u 
*J=  F(  10*6.2832 
W'.  =  »a» 

DO  1-.-  I=]*NRR 
DO  ! 5-  J=1 ,NCB 
1 39  At  I  *  J  )  =  u  . 
l'J,  CONTINUE 

DU  2-'-  J  =  1  »  K 

I F  ( MM ( J ) • LG . 1 )  GO  T u  111 
IF  <MM( J).Eu.3>  GO  TO  113 
IF  (v.V.t  JJ.F0.4)  GO  TO  114 
GO  TO  115 
111  0<J)=C<J)** 

GO  TO  115 

113  JEO=^( J) /Lut J) 

0  u  =  1  »+»<u',uEG**2 

L  C  J)  =E  (  J  )+to(  J)  **5*ufc.w**2/w.6 
•  J(  J )  =u  (  J  )  * W/GxS 
GO  TO  H5 

114  FOG  =  F'>C  J)  **?+J(  J  )  «*2 
ft J)=F( J)+FG( J)/FOG 

•J(  J)  =  F-  (j)  *<2*U(u:  /Et.o 

115  AR  =  .7o54*i)i  C  J)*i,I  ( J) 

AF I J l =AR  *E ( J ) 
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AU  (  J  )  =  AR  K  -J  (  J  ) 

DP  (  J )  =  oUi<  1  (  K  (  J )  /i\u  (  j  )  ) 

V/>P  (  J  )  =  /.'/-P  (  J  ) 
rtSP'.rdSPI  J)**? 

L  :>  (  J  )  =  L  (  J  )  *  t  l  •  +  X  (  J  )  ) 

;<r  c  j  )  =o(  j  >  /  r  f  j  > 

OCj=OE( J>**2 

COP  (  J  )  =  j  .  9  *CL>  (  J  )  *  CAP  (  J  ) 

CDK  (  J  )  =?.S284*D1  <  J  )  *L>CI<(  J)*Cl>I  /AE.1  J  ) 
o£a  =  vjE«  J)  *CDl<(  J)*.v 
GA=  a’SPS-QF  W 
GAS=GA*GA 

GB  =  CDP  (  J  )  fcA+ur(  J  >*  a’SPo 

GB3=GP*G3 

GC  =  ->uRT  ( Goj+Glli' ) 

GO =2 •  * ( 1 .+UlS) 

ALP(  J)=oCf.%T  (  (  GA  +  GC  )  /Cl> ) 

BF1  (  J)  =-ov-KI  {  ABO  (  -GA+GC  )/0U) 
ALPL=ALP( J)*LS( J) 

Bf:  U-  =  BE  T  ( J  )  *  LG  (  J ) 

.j  I  N!  1=  (  TXP  <i‘  I*  T  L  )  **(:  XP  (  -  BR  TL  )  )  /  2  • 

COGH=(  EXPdiFTL)  +(•  XP  (  “ iif  T L  )  )  /2. 

j.j  (  J  )  =GI  N(  AL.PL  )  »o  INH 

j_<  J)=GlidALPL)<C^oH 

CG(  J  )  =  COM  ALPL  )*olMI 

CC(  J)=CU.O(  ALPL)*CGb<i 

AFA( J)=AF< J)*ALP( J) 

AEP< J)=AF< J)»PPl ( J) 

AOAl  J)=A0U)»AlP<  J) 

4GB ( J ) = AO l J ) *  Br I ( J ) 

CONTINUE 

DO  o  I  =  I  *  JC 

J  =  4*  1-3  t  JJ=  J+l 

A(J»J)=GC<I)  +  A(~»J*oj)=.>C( 1 ) 

FONT  IN'JF 

no  a  v  i  =  l  *  < 

J=4*l-2  1  JJ=J-1 

A(  J»JJ) -Co  (  I  )  *  AUJ*J)---fo(I) 

COM 7  I  YJF 

IF  1  O •  L l* •  -  )  CU  10  ’vj3u 
DO  ^  vj  1=1 

Jr/**  I  -  3  i,  JJrJ+2 

A  (  J*  JJ)  =CC  {  I  )  -  A( J+1*JJ+1 )=CC<  I 

CONTINUE 

DO  6—  1  =  1  »KD 

JJ  =  4*d  -  J=JJ-3 

A  I  J  »  J  J )  =o-»  (  I  )  j>  A I  J  +  i  )  =-oo(  I  ) 

CONI  I  NuL' 

DO  7;-  I  =  1  »<D 

J  =  4  *  I  +  1  i>  JJ  =  J~2 

A(J«JJ)=-1,  *  A ( J+l *OJ+l ) =“1 • 

CON  I I  NUF 
DO  fi<--  1  =  1  *KD 
J  =  4*I~1  i>  J  K  =  J  +  1 

10=1+1 

A ( J*J~2) =Gubl 1  ) -AT A  >  i  ) 


JJ=J-1 


A  <  JJ  »  J  )  -~<  ~>{  I  } 


JJrJ+2 


A ( J+1»JJ+1 }=CC  (  I  ) 


J= JJ-3 

o> A 1  J+l  *  jJ-i  )  =  _o  j  (  I  j 


JJ= J-2 

-t  A  (  J+  3  *00+]  )  =~1  , 


J  A  =  J  +  1 


^^I>- .■  thj  s23^S3S5S5®$?#2^ 


A(  JtJ-l)=  Ann  1  )  f  A  u  A  (  !  ) 

A  (  J  ♦  J  )  =-*Vi<  I  )  **o 
A(  J>  jh  )  =-rnP(  i  j  * w 

A(J»J+2)=(ArA<I+l  )-A-J!’-(  m  )  )  *f  C  (  IM  )  »  ( AMJ(  I  U  )+AuA(  I  +  1  )  >*SS<  I  +1  ) 
A(  j.  J  +  3)  =  (  AFA  (  I  +  1  )-Awi5(  I  ♦  1  }  )  "o.'j(  1  +  1  )-(AEU(  1+1  )+AUA(  I  +  1 )  )*CC(  I  +1  ) 
A< J+1.J“2)=-A( J.J-1) 

A ( J+l * J-l ) -A< JtJ-2 ) 

A ( J+ 1 »  J  )=-A(J*J+l) 

A(J+i«J+l)=A(J*J) 


A<  J+1=J+2)  =  (AEH  (  I  u  )  +  AwA  (  Iu))’*CC(I‘-;-»  ;Awb(  Iu)-AhA(  1U)  )*bb(  Iu) 
A  (  J-*- 1  *  J  +  3  )  =  (  A'  ,Ji  I  U )  fA'^AI  I  u )  )  *bb  (  Iu ; - ( Aub (  I  u )  -AEA  (  Iu )  )*CC(  IOJ 
?.'  0  CONTINUE 

IF  (KDD.LI'.O)  GO  TO  103u 
00  ri  >o  1  =  1  tKO:) 

J  =  A  *  I  'fc  1-J=I  +  1 
JK= J-l 

\  (  Jk  »  J  f  3  )  =  (  Aft;  (  Iu  )  +  AuA(  1^  )  IHot  !■->)  +  {  A^d(  1  u  )  ~At  A  ( lu) ) *bC ( lu) 

A  (  Jr-.  »  J  +  4  )  =  (  AEH  (  Iu)+AuA(  Iu )  )  *bC  (  lu)-(  Awb(  IU)~AEA(  Iu)  )*CM  Iu) 
M  J» J+3)  =  (  -  At  A  (  10  )  +AuB<  I'-')  )Hj(  1  u  )  -  (  Ahd  (  Iu)  +AwA  (  1  u)  )  *oC  (  Iu  ) 

A ( J  »  J+4 ) - (  “At  A  (  Iu  )+Aufj(  Iu  )  )  ->bC  <  lO)  +  (  A£B(  Io)+AwA(  Iu)  )*Cb(  Iu) 
CONTINUE 
lu3u  CONTINUE 

00  1»50  I=1*NRB 

d ( I ) =u. 

1  •- 5v  CONTINUE 

n  1 1  j  =  i  • 

N  =  MRB 
M=1 

I  JUE  =  38 
JuUE  =  NCiJ 

CALL  I  .MVi<  (  A  *  ix  *  Li  *  *'i»0E  U  <m-i»  I  ul  ic»JjI  Ah  ) 

K  4  =  4  *  K  i>  A.  3  =  K4—  1 

*'J  (  <  4  )  =  0  •  i<  ti  (  K  3  )  =  u  • 

r  PRINT  5'»  ( B ( I ) » I =1 »A4 ) 

5  .  FORMAT  ( l?Elu.?) 

N  =  1 

L 1 ( 1 ) =0, 

no  1 1 uo  i  =  l » k 
KI  (  I  )=fi 
JD  =  - 

J  ,<  =  4  *  1  2  JJ  =  Jk-\ 

JI=JA~2  i  JH  =  J.<-3 

C  1  =  »'  I  J I  )  *0C  I  I  )  +'H  J<  )  *CC  *  1  )+0<  JH)  *Cb(  I  )-!>(  JJ  )*bb(  I  ) 

C?=!:<  JJ)*CCU  )+»(  J*)*oo(  1  )+:i(  JH)*SCl  I  )“B(  JI  )*CS(  I  ) 
uA(N) =S0PT (Cl =C1+C2»C2) 

VOPFC  (  N  )  =  .vb*UA  (  h  )  v*?  *  UAS  (  I  •'• ) 

Vv/AMN)  =7VAK(N)+\/bRf:C  *  M  ) 

IF  (NiC.K.Eu.0)  GO  IU  18  JO 
_  jPLC  (  N  )  =OA  (  H  )  **2*'jAb<  I  A  ) 
iJVAR  (  N  )  =UVAk  (  N  )  +U5REC  (  N  ) 

AOG=AFA( I ) *0S( I ) 

A C C  =  A [.  A  (  I  )  *rr  <  I  ) 

;.  -,r  =  AEA(  I  )  «  bC  (  I  ) 

A  c  0  =  A  L  A  (  1  )  «  r  *>  (  I  ) 

».jo=Arn(  i )  «j.>i  i ) 

:)  jC  =  Afc  :i  (  I  >  *  bC  I  I  ) 
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‘^’iTS  £Jip*3h0g3ti*e  »s 


i'.OC  =  At!U  I  )  «C  <  (  I  ) 
i*.C->  =  AhH(  I  )  *C.‘>  (  I  ) 

CArh  (  JI)  *  (  ASS-MOO  )  +iU  Jn  )  *  ( i)  ,(  4  ACS  )  -flit  JH  )*  (Mob-4  ACC  )+H(  JJ  )  *  (  ii(S-Aj>C  ) 

r  3-m (  ji  j  ••  (  ac  r  -mss » +i-. (  jk ) *  ( hc  s-aso  hm(  j.'d*  (  hoc-ass j-iu  jj  >  *  t  bso+acs  ) 
^s-r/.-oE(  1 )  *o 
0  6=0  14  or  (  I  )  »0A 

BB-ILO(N)=oul<l  (03*0  3  +  0/4*04  )  /Ah  (  I  ) 

TSHi  C  (  N  )  -U>GlLo(i4  )*  “  2*UAo(  I  u  ) 

EVAk  (  14  )  =LVAK{  14 )  4-hSPCC  ( -4) 

G  I  Gi'iA  (  N)  =  oUl<T  (05*05+06*06  ) 

T  SPFC  (  N )  =  •>  I  GHA  ( >'i )  **2*UAS  (  1 ) 

T  V  Ai<  ( N  )  =T  \CAK  ( i\  )  +T SPFC  ( N  ) 

1  M?  (  N)  =T!-12  ( .\)  +TSHEC  (  N)  *.vS 
T  M A  (  N )  =  T-V,  (  M )  +TSPEC  ( N )  *  wo*  WS 

C  hkikT  lil*ALK(  I  )  »BET(  I  )  »X.C3»CH»C>*CO»At(  I  )  »wE(  I  )  » i> I G*'*A ( ^ )  »-4 
180o  IF  (MC(  I  I.r.O.l  )  GO  TO  1210 
XINC=XINCN 
GO  TO  12?- 
1?10  XIi4C  =  XINC-> 

1220  X  I  NO  S  =  X I  NO* ( 1 «  +2 (  I  )  ) 

IF  (LI  I  ).LF,.XINC  )  GO  TO  12o- 
J  0=  I  F  I  X  (  L  (  I  )  /  X  I  i\C  ) 

DO  lib-  J= 1 » JO 
N  =  N+ 1 

LT  (i'i)=LT  ( i4- 1 )  +X  I NCS 
X  =  l.S(  I  )-Xlf4CS*FLuAl  (  J) 

ALPX  =  ALP(  I  ) «  X 
PETX  =  P.FT  (  I  )  »X 

EXH=FXP( BLTX \  i  EXBN=rXP(-UETX) 

GINhX=(FX-j-FXBN)/2. 

COSHX=  (  EXd  +  EXBiN )  /  2  • 

SoX=.3IN(  ALPX)*SI«4HX 
JCX  =  oIiM(  ALPX  )  *Cw£>MX 
C  GX  =  C>Jo  (  ALPX)  *b  I  i4HX 
CCX=C0S( ALPX)«COSHX 

C  1  =  R  (  JI  )  *->CX  +  B ( JK ) *  CCA  +B(oH)*CSX  -B(JJ)*ioA 

C?  =  vMJJ)*oCX  +B(JM*SbX  4- (  — H)  *  SC  X  -b(Jl)*CoX 
OA(N)=oOI<I  (01*01+02*02) 

V -PEC  (  N)  =..o*UA  (  M )  **2*o AS  1 1  .s  ) 

\/VAi<  ( 14 )  =\JJ AK  ( •■< )  +  VGPEC  ( 14 ) 

IF'  ( i'JK.iC  •  E'J  •  0  )  0—  To  1 1  5-- 

'JoPEC  (  N )  =-A  ( >4 )  «*2*oA3(  IK) 

Uv/A*  (  N  )  =UVAk  ( f. )  +oGPEO  (  N  ) 

ASo=ALP(  I  )  *oGX 
ACC=ALP( I )*COX 
ASC  =  ALP(  I  )*SOX 
ACS  =  ALP(  I  )  *CoX 
3JS  =  BET (  I  )  *5oa 
B sC  =  BET (  I  )*GCX 
JCC=BET ( I ) «CCX 

rcs=9et(  i  )  *csx 

C 3  =  15  C  JI  )  *(  AOC+PSo)  +  B  (  J-s)  “  (.iCG-AiO  )+b(  JH)*  (  BCC-Ao5 ) ~B ( JJ) *  (  QsC+AC-  ) 
C'4  =  o(  JI  )«(  ASS-UCC  )+o(  JX  )*(.joC  +  AoS)+d(  Jn)*(obs  +  ALC)+o(  JJ)*(uCo~Abo  ) 
C  5  =C  A-OF (  i  ) «03 
C  6=0  3  +  QE  (  I  «  *C'* 

EBol  1.0  (N  )  =  SU«  I  (  03*03+04*04  > 


H 
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I  >Pfc  C  ( ft )  =  Lt>‘>  I  LU  ( 14  )  *  a  2*U/*o  i  In) 

FVAMNj^EVAKtft)  USPI.C  (ft) 

')  I  C,A\  (  N  )  =  '->  tJ  C  T  (rr>«Cf>  +  C'A*0.»)*Ar(  I  ) 

r si»r  c<  N)  =:. igma  ( im)  **?hjas (  i  n  ) 

rVAK(.'J)  =TVAK(ft)  +1  SPFC  (ft) 
r  A?  (  im  )  =  T A  2  (  ft  )  +  r.V'EC  (  M  «  Wo 
I  ,«./*  ( IM )  =  T  .-U*  (  n  )  +  TOP  EC  <  '* )  *  fi  O *  Wo 

PRINT  l.Jl*ALP(  I  )  »UE  f  (  l  )  »X»Cb»C‘*»C:>*CO*At(  I  )  ♦«£<  l  )  »i>IG''«A("0  ,<• 

Ilf)'-  CONTINUE 
1  2fi>  ft-N+1 

K  F  (  I  )  =  \* 

ft JO=N-JD- 1 

L  T  (  N  )  •=  L  3  (  I  )+LT(NJD) 

UA  (N)  =  SORT  <  6<  JA)  *3(  Jn)  -s o<  JJ)*o<  JJ  )  ) 

VoPFC  (  N )  =*o*«JA  (  ft )  **2*lAo<  In  ) 

Vv/AKd'*)  =VVAK(u)+V5PFCCm) 

IF  (  imKA  .  f  w  •  -j  )  GU  To  18  lu 
_,.>P  EC  (  ft  )  =oA  (  ft  )**2*lJ AS ( I n ) 

OVA''  (  ft )  =UVAK  (  im  )  +OoPEC  <  ft  ) 

C3=f3  (  JI  )  *ALP(  I  )+R(  JH)*BE(  (  i  ) 

C4  =  -fU  J  I  )  *UCT  (  I  )+ij(  JH)*ALP(  1  > 

C  -=C4-0E ( 1 )*C3 
C6=C3+GE ( 1 )*C4 
Eiii>  I  LO  (  ft  )  -SON  I  (  C3*C3+C4*C4  ) 

C->fttC(ft)=LUGlLO(ii)**2’fJAS(  In) 

£\/Ai<  (  ft )  =  t“  v/AK  {  ft )  +  ESPF.C  (  a  ) 
bIGMA(N)=bUl<T(CJ)*C5+C6*C6)*AI  (  I  ) 

1  3Pf  C  (  ft )  =I>  I  Cii'iA  (  ft  )  **  ?*UAS  (  IN) 

T  VAi<  (  ft  )  =  T  V A  ft  (  A )  + 1  SPEC  (  ft  ) 

I  .-’p  (  ft  )  =TP.2  (  ft  )  +  T  SPFC  (  ft)*«MS 
I  Au  (  ft  )  =  T.'Ui  (  ft )  +T  SPEC  <  <m  )  *WS*  Wo 

PR  I  ft  T  1  b  1  *  ALP (  !  )»8ET  (  I  )  »XtC3»C,*»C3»C6*AC<  1  )  ♦'■'£(  1  )  ♦blCwAlft)  »"• 

1811'  ft  =  ft+l 

L  T  (  ft  )  =L  T  (>'•"!  HLC  I  ) 

H--*  CONTINUE 

IF  (IiKK.Fm.u)  GJ  TO  182^ 

4 ■  ■  1 1  /  NH=  M— 1 

PK  IN  r  6-  *  F  (  IK) 

60  Foi'..,AT  (  1  r  i  »/////////k  iHIo  lo  I  -)L  RtoUi-  1  b  v<  1 1  H  D  A‘"‘P  I  ‘"G  AT  F*m£  >» 
loc-CYr*  »Fc.6/////9X»',oA”  » 1  OX  »  ’‘Uo^LC*  »  10X  »  *Eoo  I  »9X  » ■•‘LSi'LC’'  *  1UX  *  * 

loIG'-'A*  »luX,«T  >PEC*  » 12X  »  *0  I  *,  1 }A**I*  ) 
p>'  IAT70*  (uA(  I  )  *USPCC(  I  )  »Cc>oI«-o(  I  )  ,  EbPLCI  1  )  »o  ’.'onAt  I  )  *TSPEC(  1  )  *lI  ( 
1I)*I»  I  -  1 » N'H ) 

70  FORMAT (7E1 3.4* I lo) 

18?  •  -< E ^  F  C  r.  ) 

PRINT  133*  (UA  ( J  )  *  Vo^EC  <-• )  »^vA>'(  j  )  .uT  (o  )  ,  o  =  l»  -'£n) 

its  fqp.-'.at (4E14.4) 

.  CONTINUE 

PK I  ft  I  13?*  (v/VA<(0)»  LT(o)»  J=1»n£^) 

]  3?  FORMAT  ( 2. F  14.4) 

PRINT  1 33 » < KJ ( J) »<E ( J ) *  J=1*n) 

188  FORMAT  ( ? I  1 o ) 

IF  riiCK.E-.U)  GO  TO  183 
■"»o  ? - ■  1  r,  i  =  i*:ih 

I  .-iw4=  T  VAN  (  I  )*T  -’.4  C  I  ) 
r  ;•)??  =  tm?  { 1  j  *  t y>2  ( i ) 


U"' 


4'  ■  !  '/ 
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I  I  -  (  I  I'-'l  “  I  Pi.V  )/  1  MO  4 
1  l.l*.  (  I  )  -  VJl<  MIL-  H-> J 
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